c Information Technology Laboratory, National Institute of 
c Standards and Technology, 325 Broadway, Mail Stop 891.01, 
c Boulder, CO 80305-3328. email: zydrunas.gimbutas@nist.gov. 
c Contributions by staff of NIST, an agency of the U.S. Government, 
c are not subject to copyright within the United States.

c NIST-developed software is provided by NIST as a public service. 
c You may use, copy, and distribute copies of the software in any medium, 
c provided that you keep intact this entire notice. You may improve, 
c modify, and create derivative works of the software or any portion of 
c the software, and you may copy and distribute such 
c modifications or works. Modified works should carry a notice stating 
c that you changed the software and should note the date and nature 
c of any such change. Please explicitly acknowledge the National 
c Institute of Standards and Technology as the source of the software. 

c NIST-developed software is expressly provided "AS IS." 
c NIST MAKES NO WARRANTY OF ANY KIND, EXPRESS, IMPLIED, IN FACT, OR 
c ARISING BY OPERATION OF LAW, INCLUDING, WITHOUT LIMITATION, THE 
c IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE, 
c NON-INFRINGEMENT, AND DATA ACCURACY. NIST NEITHER REPRESENTS NOR 
c WARRANTS THAT THE OPERATION OF THE SOFTWARE WILL BE UNINTERRUPTED 
c OR ERROR-FREE, OR THAT ANY DEFECTS WILL BE CORRECTED. NIST DOES NOT
c WARRANT OR MAKE ANY REPRESENTATIONS REGARDING THE USE OF THE 
c SOFTWARE OR THE RESULTS THEREOF, INCLUDING BUT NOT LIMITED TO THE 
c CORRECTNESS, ACCURACY, RELIABILITY, OR USEFULNESS OF THE SOFTWARE.
c
c You are solely responsible for determining the appropriateness of 
c using and distributing the software and you assume all risks 
c associated with its use, including but not limited to the risks 
c and costs of program errors, compliance with applicable laws, damage 
c to or loss of data, programs or equipment, and the unavailability 
c or interruption of operation. This software is not intended to be 
c used in any situation where a failure could cause risk of injury 
c or damage to property. The software developed by NIST employees is 
c not subject to copyright protection within the United States.


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c

c     this file contains plane wave quadrature info
c      for sepration by 1
c    
      subroutine vwts(x,w,n)
	  implicit real*8 (a-h,o-z)
      dimension w(n),x(n)
c
c
c This routine returns a set of Gaussian nodes and weights for
c integrating the functions J0(r*x)*exp(z*x)dx over the range x=0 to
c x=infinity.  
c They work only for values of z within the range [1,4] and r within
c the range [0,4*sqrt(2)].
c
c
c Input arguments:
c	n - number of weights and nodes in the quadrature.  This must
c	  be an integer in the range [2,39].
c
c Output arguments:
c	w - weights
c	x - nodes
c
c
c Approximate accuracies of the quadrature:
c
c	N	maximum error
c
c        2     0.15318D+00
c        3     0.76505D-01
c        4     0.32149D-01
c        5     0.15630D-01
c        6     0.75110D-02
c        7     0.35030D-02
c        8     0.16243D-02
c        9     0.72230D-03
c       10     0.33074D-03
c       11     0.15035D-03
c       12     0.70952D-04
c       13     0.31751D-04
c       14     0.14589D-04
c       15     0.64300D-05
c       16     0.29477D-05
c       17     0.13222D-05
c       18     0.61488D-06
c       19     0.27435D-06
c       20     0.12534D-06
c       21     0.55324D-07
c       22     0.25257D-07
c       23     0.11293D-07
c       24     0.52063D-08
c       25     0.23256D-08
c       26     0.10580D-08
c       27     0.46835D-09
c       28     0.21286D-09
c       29     0.95164D-10
c       30     0.43599D-10
c       31     0.19516D-10
c       32     0.88491D-11
c       33     0.39313D-11
c       34     0.17821D-11
c       35     0.79603D-12
c       36     0.36460D-12
c       37     0.16331D-12
c       38     0.73497D-13
c       39     0.31530D-13
c
c
	if(n.lt.2 .or. n.gt.39) stop 'n out of bounds in vwts()'

	if(n.eq.  2) then
	  x(1)=0.31283512406199648347993047536874656D+00
	  x(2)=0.10894076617764414383060511681833304D+01
	  w(1)=0.76974593066577023936503110235207714D+00
	  w(2)=0.84374978799072331003827684980933554D+00
	endif
	if(n.eq.  3) then
	  x(1)=0.23817561600133754784458517406164901D+00
	  x(2)=0.90437226188379571389930333680240437D+00
	  x(3)=0.17511910669445907018371144658885896D+01
	  w(1)=0.55695724078723707339833026708220132D+00
	  w(2)=0.79531004550917117068564721193979494D+00
	  w(3)=0.93703416373547598894333532371092588D+00
	endif
	if(n.eq.  4) then
	  x(1)=0.21626008006106806069723802465887275D+00
	  x(2)=0.85186788888114162165976495089125820D+00
	  x(3)=0.16705322595388429895990611839806661D+01
	  x(4)=0.25585029089804702806532077374868095D+01
	  w(1)=0.50471599046036208502385989049798809D+00
	  w(2)=0.76071728649476022532383012730861083D+00
	  w(3)=0.87370748029129807754600278713041916D+00
	  w(4)=0.93475133313404745738495194018469192D+00
	endif
	if(n.eq.  5) then
	  x(1)=0.17368350550918876251671463251113892D+00
	  x(2)=0.71973899919337358532800408283947036D+00
	  x(3)=0.14665863986721390244127860569278710D+01
	  x(4)=0.23283854455252543935728226642822847D+01
	  x(5)=0.32597538486988750783268642408074811D+01
	  w(1)=0.40839334293432694167336194368544966D+00
	  w(2)=0.67016366313130026899358426817343570D+00
	  w(3)=0.81673864902764081730168754802434705D+00
	  w(4)=0.91106425433364368604571836840477772D+00
	  w(5)=0.97880305318787497181887147235102020D+00
	endif
	if(n.eq.  6) then
	  x(1)=0.14624625076566644121811577861080877D+00
	  x(2)=0.63832297842729379677706447182572447D+00
	  x(3)=0.13388327339170913443666677267174236D+01
	  x(4)=0.21582290084516078110254966304637492D+01
	  x(5)=0.30493033990054105863976019463734701D+01
	  x(6)=0.39906171458431614063044889917364344D+01
	  w(1)=0.34993298510539821988274411523889285D+00
	  w(2)=0.61733000302790996816071356079191901D+00
	  w(3)=0.77195750484413483416545886939275078D+00
	  w(4)=0.86303057957197659710857351456070319D+00
	  w(5)=0.92540796896975519025829726160736755D+00
	  w(6)=0.98300735313837050988183818844845518D+00
	endif
	if(n.eq.  7) then
	  x(1)=0.12110405743821943402860341620908002D+00
	  x(2)=0.56002252219580128667075769044458866D+00
	  x(3)=0.12107447642184836489320787222823128D+01
	  x(4)=0.19939530286622262345019862550543621D+01
	  x(5)=0.28652677773020815976678932202048600D+01
	  x(6)=0.37987093416333359563452631846303120D+01
	  x(7)=0.47769841417427238994264371285680681D+01
	  w(1)=0.29698283202048353368596167456416879D+00
	  w(2)=0.56247171144261287611954003295977600D+00
	  w(3)=0.72718448482782827024806238114251755D+00
	  w(4)=0.83360054399649463618970912648364902D+00
	  w(5)=0.90833570000297958468848946722573601D+00
	  w(6)=0.96486190247307246359298460447462276D+00
	  w(7)=0.10169919877814406206084640871267766D+01
	endif
	if(n.eq.  8) then
	  x(1)=0.10934746769000448785025270126425312D+00
	  x(2)=0.51769741015341219636525238456670195D+00
	  x(3)=0.11330659161119163158559786097612232D+01
	  x(4)=0.18813501511074042227278368955012411D+01
	  x(5)=0.27178540960120538549915636394871399D+01
	  x(6)=0.36165027490744909854925026593264192D+01
	  x(7)=0.45627105330382127235111511254217476D+01
	  x(8)=0.55490088534852830548516067210584879D+01
	  w(1)=0.27107502662773547141839003415952902D+00
	  w(2)=0.52769158843946417647430280339904130D+00
	  w(3)=0.69151504413879483212923560131457634D+00
	  w(4)=0.79834400406452044052230121451430023D+00
	  w(5)=0.87164160121353972954238997772336006D+00
	  w(6)=0.92643839116924142640385753111331724D+00
	  w(7)=0.97294622259483076742725415897439234D+00
	  w(8)=0.10241386584468554765692260843934491D+01
	endif
	if(n.eq.  9) then
	  x(1)=0.99273996739714473469540223504736787D-01
	  x(2)=0.47725674637049431137114652301534079D+00
	  x(3)=0.10553366138218296388373573790886439D+01
	  x(4)=0.17675934335400844688024335482623428D+01
	  x(5)=0.25734262935147067530294862081063911D+01
	  x(6)=0.34482433920158257478760788217186928D+01
	  x(7)=0.43768098355472631055818055756390095D+01
	  x(8)=0.53489575720546005399569367000367492D+01
	  x(9)=0.63576578531337464283978988532908261D+01
	  w(1)=0.24776441819008371281185532097879332D+00
	  w(2)=0.49188566500464336872511239562300034D+00
	  w(3)=0.65378749137677805158830324216978624D+00
	  w(4)=0.76433038408784093054038066838984378D+00
	  w(5)=0.84376180565628111640563702167128213D+00
	  w(6)=0.90445883985098263213586733400006779D+00
	  w(7)=0.95378613136833456653818075210438110D+00
	  w(8)=0.99670261613218547047665651916759089D+00
	  w(9)=0.10429422730252668749528766056755558D+01
	endif
	if(n.eq. 10) then
	  x(1)=0.92181119494357380284021985517028952D-01
	  x(2)=0.44689957388438356078452784458932001D+00
	  x(3)=0.99485853124047674000962615536991507D+00
	  x(4)=0.16757226616017673048730785012594424D+01
	  x(5)=0.24510992450330819636405976780224591D+01
	  x(6)=0.32961838574354240805064364394638687D+01
	  x(7)=0.41949178113509839604944318125490099D+01
	  x(8)=0.51371729038633668906754792260471731D+01
	  x(9)=0.61167060097997145717840794532094151D+01
	  x(10)=0.71305602574974669138896388176362962D+01
	  w(1)=0.23083073738461035162572443368844688D+00
	  w(2)=0.46385202776449724337837210441648494D+00
	  w(3)=0.62235404858126219806990775396116078D+00
	  w(4)=0.73321467536730178249371192578109913D+00
	  w(5)=0.81362306975392817420811297779437155D+00
	  w(6)=0.87427212844441348771340472012525424D+00
	  w(7)=0.92246427999726521296253167747636326D+00
	  w(8)=0.96346501535401229965316360903671011D+00
	  w(9)=0.10015863043245734065322949390974827D+01
	  w(10)=0.10469665614659515728845917692524381D+01
	endif
	if(n.eq. 11) then
	  x(1)=0.83947754993541900958575752156320959D-01
	  x(2)=0.41116156884118149061890790108009242D+00
	  x(3)=0.92412911280064913288612160613411106D+00
	  x(4)=0.15697112352606923391107329734950326D+01
	  x(5)=0.23129015267322943216754538298118860D+01
	  x(6)=0.31300036645182620631544523348566145D+01
	  x(7)=0.40050323813231880620833180728368461D+01
	  x(8)=0.49273172978328521054436350823380053D+01
	  x(9)=0.58896866519437089948496577562764287D+01
	  x(10)=0.68870955302308018985968374181538820D+01
	  x(11)=0.79164395304058903235500110895372927D+01
	  w(1)=0.21096619809745237317244459518406074D+00
	  w(2)=0.43092534827878592462724327560863458D+00
	  w(3)=0.58635954971056480911073549577849917D+00
	  w(4)=0.69913461337603133305407254738383926D+00
	  w(5)=0.78339916419475752462631135131232440D+00
	  w(6)=0.84830422968296659291809191927313805D+00
	  w(7)=0.90031891956613319116797811147989705D+00
	  w(8)=0.94389939178410831033261274569667876D+00
	  w(9)=0.98221287456430605988799698025104590D+00
	  w(10)=0.10179988215575410670510336785810068D+01
	  w(11)=0.10599575740381281807600544198066927D+01
	endif
	if(n.eq. 12) then
	  x(1)=0.76855224482364675742118720336293336D-01
	  x(2)=0.38022716855965127402683378932124469D+00
	  x(3)=0.86295016672459190498756242959643714D+00
	  x(4)=0.14774065742425332992127096076728776D+01
	  x(5)=0.21905930725126028235649755515623838D+01
	  x(6)=0.29791885550546846062047734449151903D+01
	  x(7)=0.38268052131682357419606432813452557D+01
	  x(8)=0.47221812142851433691248530521988869D+01
	  x(9)=0.56578288522785102898637887847144157D+01
	  x(10)=0.66290084039626417933277480187825859D+01
	  x(11)=0.76329115194492631601974608201999217D+01
	  x(12)=0.86692585676959215135184422251768410D+01
	  w(1)=0.19378032292424979265987872167897876D+00
	  w(2)=0.40247808945013635595344680950802285D+00
	  w(3)=0.55512328548655365256792038053390570D+00
	  w(4)=0.66840122968153037330552024286589585D+00
	  w(5)=0.75414462244054158990280711805098690D+00
	  w(6)=0.82039053613530976694079299704753794D+00
	  w(7)=0.87310177781587317014100335654802620D+00
	  w(8)=0.91671095974371530346047620696481317D+00
	  w(9)=0.95457288752598756786937883589416742D+00
	  w(10)=0.98937097491594594433905740515911020D+00
	  w(11)=0.10238743680564139637567677709739655D+01
	  w(12)=0.10678249338234333531971742559107952D+01
	endif
	if(n.eq. 13) then
	  x(1)=0.70687308611330559870999934446444968D-01
	  x(2)=0.35322227031870156999104892747709528D+00
	  x(3)=0.80956306466728034099133992640418001D+00
	  x(4)=0.13969368075858923372578601629356854D+01
	  x(5)=0.20843093180753342608113598544150591D+01
	  x(6)=0.28492369286741880962665618426399305D+01
	  x(7)=0.36755557385319534979828404175350443D+01
	  x(8)=0.45518102122361847605702678265515715D+01
	  x(9)=0.54700092351619344910318432084750384D+01
	  x(10)=0.64246720832453547700424678623676300D+01
	  x(11)=0.74121271804332247157276469806674868D+01
	  x(12)=0.84299356852023183961364338756538928D+01
	  x(13)=0.94776462234869871537057406385429204D+01
	  w(1)=0.17879600056265629315355170092516346D+00
	  w(2)=0.37763057118432419301257141341920942D+00
	  w(3)=0.52791576609439927469935582848847844D+00
	  w(4)=0.64173473833746008665457338793203235D+00
	  w(5)=0.72931281728819530307816876302240416D+00
	  w(6)=0.79789045633520616629397181895910762D+00
	  w(7)=0.85289957189509402191163189854705706D+00
	  w(8)=0.89839690889403756024478298058966175D+00
	  w(9)=0.93736682912526381183937473906553350D+00
	  w(10)=0.97206402733131802573751656382228248D+00
	  w(11)=0.10043560920792620194674782396759838D+01
	  w(12)=0.10363376029237720832298919049208052D+01
	  w(13)=0.10772868663681449685753932499210350D+01
	endif
	if(n.eq. 14) then
	  x(1)=0.66213023807563053946267928040469997D-01
	  x(2)=0.33302518940702452443503034373861738D+00
	  x(3)=0.76814378543799699539107450618757866D+00
	  x(4)=0.13320738699370830993018444132758304D+01
	  x(5)=0.19952905542395280580336702769272961D+01
	  x(6)=0.27361894979680170258973248564871028D+01
	  x(7)=0.35389332083102940096352995169581845D+01
	  x(8)=0.43920171385366097993596667947713286D+01
	  x(9)=0.52872133930334470619527564849704504D+01
	  x(10)=0.62187865433268365578101111168507487D+01
	  x(11)=0.71829341316549317753015202470123768D+01
	  x(12)=0.81773810009317156044517105328850448D+01
	  x(13)=0.92010945943890938281128910603001714D+01
	  x(14)=0.10255597561090382896509254351258278D+02
	  w(1)=0.16782217858320716552533724552631611D+00
	  w(2)=0.35834278682050713094042748707579449D+00
	  w(3)=0.50524782135037205055283493493334390D+00
	  w(4)=0.61774741033421598501718108309432864D+00
	  w(5)=0.70512340043420818069819233642192557D+00
	  w(6)=0.77406410364509525834364467300474644D+00
	  w(7)=0.82953851602330785919292566177318804D+00
	  w(8)=0.87530868438704301226493953436147422D+00
	  w(9)=0.91423294730365067817245972037198953D+00
	  w(10)=0.94853481721964949535674804792506620D+00
	  w(11)=0.98004492543042398544628213130636141D+00
	  w(12)=0.10104582738217686710413545370101929D+01
	  w(13)=0.10420257056214121682558015891117975D+01
	  w(14)=0.10848830289657354075671946702641435D+01
	endif
	if(n.eq. 15) then
	  x(1)=0.62494006208776084476674839152110508D-01
	  x(2)=0.31580001757481546142969364154851064D+00
	  x(3)=0.73196261825059882877297923187143169D+00
	  x(4)=0.12745096769359307575797402023454197D+01
	  x(5)=0.19156615424369782818558860526536591D+01
	  x(6)=0.26349779278743241661686624865978956D+01
	  x(7)=0.34172786238081047649473021010635421D+01
	  x(8)=0.42512905564487253329275517899077386D+01
	  x(9)=0.51287187209510305407889063644688576D+01
	  x(10)=0.60435472353063381945048604393377900D+01
	  x(11)=0.69915157252578845614721103629563004D+01
	  x(12)=0.79697454656597592759226245107129216D+01
	  x(13)=0.89764624126137402271297105471603572D+01
	  x(14)=0.10010758126582684823802082973998040D+02
	  x(15)=0.11073815811281715326686025946401060D+02
	  w(1)=0.15862394244419067423912395042862045D+00
	  w(2)=0.34142839128318125396788218495203182D+00
	  w(3)=0.48472220331313015106999841918877792D+00
	  w(4)=0.59578598263167215076663296713377349D+00
	  w(5)=0.68314771258188122793342245131498203D+00
	  w(6)=0.75297871300677499117881552592734806D+00
	  w(7)=0.80976432197984060668005668048863299D+00
	  w(8)=0.85690189548240613959251277265138924D+00
	  w(9)=0.89699396665012665064153907223953865D+00
	  w(10)=0.93204514056719467429701353466953151D+00
	  w(11)=0.96364174051181383795494639343814924D+00
	  w(12)=0.99312923342960057748740609895321541D+00
	  w(13)=0.10218026886214661708152107166824862D+01
	  w(14)=0.10514581504829776381626516013056971D+01
	  w(15)=0.10915651423913341755422834467026405D+01
	endif
	if(n.eq. 16) then
	  x(1)=0.59273387891186600506809867283664062D-01
	  x(2)=0.30052912864803243975231339391029906D+00
	  x(3)=0.69913214785198229961338256543967873D+00
	  x(4)=0.12212613306065938090227973589207977D+01
	  x(5)=0.18407022301918960494049315457232296D+01
	  x(6)=0.25380134246042680423727233574027196D+01
	  x(7)=0.32985714197429576444164922577328980D+01
	  x(8)=0.41112958629313363090318489412311465D+01
	  x(9)=0.49678340339823749971515098877716810D+01
	  x(10)=0.58619693058738864266388191026635468D+01
	  x(11)=0.67891769853120012001568284176755697D+01
	  x(12)=0.77463043528356614331187302013859153D+01
	  x(13)=0.87313537759783503844346341793425381D+01
	  x(14)=0.97433360522222471900022355839610100D+01
	  x(15)=0.10782205113039477595293647027574480D+02
	  x(16)=0.11850487340578309769512088678311557D+02
	  w(1)=0.15059862886412836435745532526198076D+00
	  w(2)=0.32605002595383431041042854303668719D+00
	  w(3)=0.46540860398333899761169618614076171D+00
	  w(4)=0.57451518283541402887237836694112048D+00
	  w(5)=0.66116138075584618416513649208354764D+00
	  w(6)=0.73104478829065544509546725748805329D+00
	  w(7)=0.78823951974329575786981649798690341D+00
	  w(8)=0.83583181691010888236803566542221233D+00
	  w(9)=0.87622783297370820942973068667924963D+00
	  w(10)=0.91132442140861513735217158682644367D+00
	  w(11)=0.94264787385567350153081633834517561D+00
	  w(12)=0.97148379898864845038986004510661587D+00
	  w(13)=0.99900986222142706960624991552322172D+00
	  w(14)=0.10264744781310688459541324846213683D+01
	  w(15)=0.10558364249820026969928221660666168D+01
	  w(16)=0.10972184518756178395193501273752190D+01
	endif
	if(n.eq. 17) then
	  x(1)=0.55990025317491039225092919195958530D-01
	  x(2)=0.28485138101968032753319448602269404D+00
	  x(3)=0.66535367065853057155777605657931417D+00
	  x(4)=0.11666790480529627416217408608645201D+01
	  x(5)=0.17644302741343080853653191297780722D+01
	  x(6)=0.24402983223637999543598198215477169D+01
	  x(7)=0.31803218099151462361362519004615024D+01
	  x(8)=0.39737171577719276882589838351123035D+01
	  x(9)=0.48121679941063382202059983683284372D+01
	  x(10)=0.56893231451148746913304421468637884D+01
	  x(11)=0.66004047944437740369494349579326808D+01
	  x(12)=0.75419049746991051463851363223511726D+01
	  x(13)=0.85113656929809931028785285889171064D+01
	  x(14)=0.95072324275912780677799673867411911D+01
	  x(15)=0.10528748096509671583476119849365205D+02
	  x(16)=0.11575870196028841263569120201282203D+02
	  x(17)=0.12650781639685202861755897174589336D+02
	  w(1)=0.14239483712193726216987954558135243D+00
	  w(2)=0.31017671029270804661237548316421453D+00
	  w(3)=0.44557516683708542348796299847890623D+00
	  w(4)=0.55303383994158716863154268139624037D+00
	  w(5)=0.63944903363522653805262052628677338D+00
	  w(6)=0.70997911214019127079666304780403152D+00
	  w(7)=0.76828253949732216998569356292136945D+00
	  w(8)=0.81713201141707436647720896871760488D+00
	  w(9)=0.85872191623337057464482313662301749D+00
	  w(10)=0.89480789582390418690494016118464060D+00
	  w(11)=0.92680189417317260236472975520882756D+00
	  w(12)=0.95586282708095549054405637434683740D+00
	  w(13)=0.98299145008229538156285798322642222D+00
	  w(14)=0.10091339538570300859987582953181118D+01
	  w(15)=0.10353177460050846381278688568272628D+01
	  w(16)=0.10631842791396344249932326420093887D+01
	  w(17)=0.11023210952108835058993463462684304D+01
	endif
	if(n.eq. 18) then
	  x(1)=0.52788527661177607475107009804560221D-01
	  x(2)=0.26949859838931256028615734976483509D+00
	  x(3)=0.63220353174689392083962502510985360D+00
	  x(4)=0.11130756427760852833586113774799742D+01
	  x(5)=0.16893949614021379623807206371566281D+01
	  x(6)=0.23437620046953044905535534780938178D+01
	  x(7)=0.30626998290780611533534738555317745D+01
	  x(8)=0.38356294126529686394633245072327554D+01
	  x(9)=0.46542473432156272750148673367220908D+01
	  x(10)=0.55120938659358147404532246582675725D+01
	  x(11)=0.64042126837727888499784967279992998D+01
	  x(12)=0.73268800190617540124549122992902994D+01
	  x(13)=0.82774009925823861522076185792684555D+01
	  x(14)=0.92539718060248947750778825138695538D+01
	  x(15)=0.10255602723746401139237605093512684D+02
	  x(16)=0.11282088297877740146191172243561596D+02
	  x(17)=0.12334067909676926788620221486780792D+02
	  x(18)=0.13414920240172401477707353478763252D+02
	  w(1)=0.13438265914335215112096477696468355D+00
	  w(2)=0.29457752727395436487256574764614925D+00
	  w(3)=0.42607819361148618897416895379137713D+00
	  w(4)=0.53189220776549905878027857397682965D+00
	  w(5)=0.61787306245538586857435348065337166D+00
	  w(6)=0.68863156078905074508611505734734237D+00
	  w(7)=0.74749099381426187260757387775811367D+00
	  w(8)=0.79699192718599998208617307682288811D+00
	  w(9)=0.83917454386997591964103548889397644D+00
	  w(10)=0.87570092283745315508980411323136650D+00
	  w(11)=0.90792943590067498593754180546966381D+00
	  w(12)=0.93698393742461816291466902839601971D+00
	  w(13)=0.96382546688788062194674921556725167D+00
	  w(14)=0.98932985769673820186653756536543369D+00
	  w(15)=0.10143828459791703888726033255807124D+01
	  w(16)=0.10400365437416452252250564924906939D+01
	  w(17)=0.10681548926956736522697610780596733D+01
	  w(18)=0.11090758097553685690428437737864442D+01
	endif
	if(n.eq. 19) then
	  x(1)=0.49949579666168537894765933060625684D-01
	  x(2)=0.25587785751621588969939580238133203D+00
	  x(3)=0.60285193295998951157343981321901083D+00
	  x(4)=0.10657983137580619992235142490244471D+01
	  x(5)=0.16235028084630112044095540113630705D+01
	  x(6)=0.22593565245709532618434423056896776D+01
	  x(7)=0.29603306461544551808628966682590544D+01
	  x(8)=0.37160847801044232241451936715748161D+01
	  x(9)=0.45183959543781595513678439601790160D+01
	  x(10)=0.53607645945685540311842487426474690D+01
	  x(11)=0.62381054190685949478734073636587709D+01
	  x(12)=0.71464949862225219945344178995583206D+01
	  x(13)=0.80829720193041367792829987592995167D+01
	  x(14)=0.90453910811754347776059148600324988D+01
	  x(15)=0.10032329160512597709953297453466803D+02
	  x(16)=0.11043040453225657415714522358030081D+02
	  x(17)=0.12077441113209841816455991647671908D+02
	  x(18)=0.13136150656231230016146582784131169D+02
	  x(19)=0.14222379224644413397982134483754635D+02
	  w(1)=0.12727479585254397598248488066019490D+00
	  w(2)=0.28074654471419330237225153723557014D+00
	  w(3)=0.40890957644635866907734111919126008D+00
	  w(4)=0.51344857052912906691943817349965684D+00
	  w(5)=0.59921282467118286430007856324664317D+00
	  w(6)=0.67033798789497855974417461766279303D+00
	  w(7)=0.72989603646689105165279443099279888D+00
	  w(8)=0.78024834204579218699393550195964053D+00
	  w(9)=0.82329643252275364062597873271442950D+00
	  w(10)=0.86059894129480551949740174677572213D+00
	  w(11)=0.89343570954146100948634057203889824D+00
	  w(12)=0.92286139944663625467313750050379895D+00
	  w(13)=0.94976012738946247182525439711753279D+00
	  w(14)=0.97490342742485036620081473301979713D+00
	  w(15)=0.99901486587368271141684772373992018D+00
	  w(16)=0.10228499282064815734116791645647027D+01
	  w(17)=0.10473249251926426417469428997719660D+01
	  w(18)=0.10741045635513459810539416139363311D+01
	  w(19)=0.11134033950443573690591847480391152D+01
	endif
	if(n.eq. 20) then
	  x(1)=0.47543727637248914263956578452052781D-01
	  x(2)=0.24422020867640242824947449662431609D+00
	  x(3)=0.57738735420750131321909748294274323D+00
	  x(4)=0.10241481132215648486294412577990443D+01
	  x(5)=0.15644872035736114845860811328748241D+01
	  x(6)=0.21824433236293536708672036183997989D+01
	  x(7)=0.28653954538190715695122889883350581D+01
	  x(8)=0.36032679528582369776756877399748191D+01
	  x(9)=0.43879813157233691356395866023376584D+01
	  x(10)=0.52130765667652978478940895001869649D+01
	  x(11)=0.60734306526893613309425745683256537D+01
	  x(12)=0.69650271185981909027873371087480336D+01
	  x(13)=0.78847726403958340313238295493647456D+01
	  x(14)=0.88303579521578470234999258536845446D+01
	  x(15)=0.98001630879826002740173862548545003D+01
	  x(16)=0.10793207610765842119349144923035055D+02
	  x(17)=0.11809146586022226443901672610081732D+02
	  x(18)=0.12848309996690652212691929889842868D+02
	  x(19)=0.13911859887638442501156532671302557D+02
	  x(20)=0.15003974288984196050478203687816858D+02
	  w(1)=0.12123489815809705694160669509074069D+00
	  w(2)=0.26876608765865678263295990291226190D+00
	  w(3)=0.39364274681026828472951706316962373D+00
	  w(4)=0.49653435789152050583439290676324163D+00
	  w(5)=0.58149909773786201405698648159159347D+00
	  w(6)=0.65232326071962210001231596834259108D+00
	  w(7)=0.71191063049174330146229294769000262D+00
	  w(8)=0.76249393422317224278117464564274997D+00
	  w(9)=0.80586118795410766058040508141857572D+00
	  w(10)=0.84347981485805434420655046778847463D+00
	  w(11)=0.87656332288607641256561464615515433D+00
	  w(12)=0.90611967404999338437932010492659174D+00
	  w(13)=0.93299649036540077684520611001062207D+00
	  w(14)=0.95792683395701161419566460608621128D+00
	  w(15)=0.98157794960527799954519423408783041D+00
	  w(16)=0.10046090612120206486679308000020683D+01
	  w(17)=0.10277523675616564435841837621410377D+01
	  w(18)=0.10519663093854307867047737090615556D+01
	  w(19)=0.10791179202989034102699861250584945D+01
	  w(20)=0.11199933499526084634823064334341325D+01
	endif
	if(n.eq. 21) then
	  x(1)=0.45571106583094517228360587068891618D-01
	  x(2)=0.23456924197771603091133840734983096D+00
	  x(3)=0.55605074355978634237374080839799717D+00
	  x(4)=0.98886226211903260896463052631588653D+00
	  x(5)=0.15140517506819852755484134831931442D+01
	  x(6)=0.21163084058946690646507704514078796D+01
	  x(7)=0.27834664049554236342487456568051130D+01
	  x(8)=0.35057749223342495881183822348248214D+01
	  x(9)=0.42753582277771142727829101204406470D+01
	  x(10)=0.50858504213438919805412297137081623D+01
	  x(11)=0.59321299588987201900636136997491121D+01
	  x(12)=0.68101102326376521389761364844162017D+01
	  x(13)=0.77165699408569325967732766002882272D+01
	  x(14)=0.86490189544857720704840176040306687D+01
	  x(15)=0.96055990356413225583764869952574372D+01
	  x(16)=0.10585019994373375595841935137286782D+02
	  x(17)=0.11586532698006379860089509747922421D+02
	  x(18)=0.12609939047089179453564611321780831D+02
	  x(19)=0.13655630902209903609900720766745508D+02
	  x(20)=0.14724711974343017928390509041491896D+02
	  x(21)=0.15821115878987425773516406479757279D+02
	  w(1)=0.11626932798637451549339516532199923D+00
	  w(2)=0.25873576304378226309310662145435344D+00
	  w(3)=0.38059684744358213137616075982805341D+00
	  w(4)=0.48185261255750905240446968491596635D+00
	  w(5)=0.56598645649837764359801894897827879D+00
	  w(6)=0.63651204484922907500532573976670392D+00
	  w(7)=0.69618766407552234998235007878975011D+00
	  w(8)=0.74712267140691357436566022443003021D+00
	  w(9)=0.79098736463667274065869605692569166D+00
	  w(10)=0.82914708606812959068577129073673859D+00
	  w(11)=0.86273418493909082016557476890739053D+00
	  w(12)=0.89269289742930946207621900612139143D+00
	  w(13)=0.91981614271009970451586923445574939D+00
	  w(14)=0.94478093124118001089328799935174175D+00
	  w(15)=0.96818477444619643890177940193098038D+00
	  w(16)=0.99058577205436648682734812609851360D+00
	  w(17)=0.10125517975529727543460012384457514D+01
	  w(18)=0.10347297124769503362529121659463271D+01
	  w(19)=0.10579732423280954201061376807047054D+01
	  w(20)=0.10839601118832196124941447124001570D+01
	  w(21)=0.11232230992403440250626545093837194D+01
	endif
	if(n.eq. 22) then
	  x(1)=0.43730571482595376187063607176241931D-01
	  x(2)=0.22546314061070582557277930391137488D+00
	  x(3)=0.53562452322468856458215213933726773D+00
	  x(4)=0.95459967559048286567247032508021221D+00
	  x(5)=0.14644503885175594870560189519892447D+01
	  x(6)=0.20505129898819740574822390044573694D+01
	  x(7)=0.27010747244419337498300137667683885D+01
	  x(8)=0.34067100069124638395123838563449681D+01
	  x(9)=0.41597567284511374197109034867025912D+01
	  x(10)=0.49539626251771764486875326838344336D+01
	  x(11)=0.57842378614529854630177396757062525D+01
	  x(12)=0.66464658921597612462051074544433504D+01
	  x(13)=0.75373509780144951619718085567001253D+01
	  x(14)=0.84542953579871529257161455461755395D+01
	  x(15)=0.93953045914506567726220964686945081D+01
	  x(16)=0.10358921378680195246602124825585634D+02
	  x(17)=0.11344189462364381526526813104283065D+02
	  x(18)=0.12350651043352140590059207170270383D+02
	  x(19)=0.13378382354647602170416575972922146D+02
	  x(20)=0.14428071458844170749102886475156993D+02
	  x(21)=0.15501228719971146929879068920854479D+02
	  x(22)=0.16602530548254346598469055606983602D+02
	  w(1)=0.11162205199242761222144793009647401D+00
	  w(2)=0.24914697542128799723215593076020014D+00
	  w(3)=0.36780722045157282495964068402827252D+00
	  w(4)=0.46713432561546514021344478351238649D+00
	  w(5)=0.55013206743428866296596879692515358D+00
	  w(6)=0.62005292437552328443928217893699184D+00
	  w(7)=0.67950829126698175652876443564309739D+00
	  w(8)=0.73049097850935695852570006536552683D+00
	  w(9)=0.77456446199081319115009591769194230D+00
	  w(10)=0.81300220581189297508473146081087179D+00
	  w(11)=0.84686488163552398056310721585759893D+00
	  w(12)=0.87704436975849753288514421001309529D+00
	  w(13)=0.90429550558952997985073807285516523D+00
	  w(14)=0.92926471016166922733958699609502219D+00
	  w(15)=0.95251860675498312147624346835073084D+00
	  w(16)=0.97457439022305569764625943207647651D+00
	  w(17)=0.99593510961900155731285622096038423D+00
	  w(18)=0.10171365736308934923926017290796153D+01
	  w(19)=0.10388201977843312118210405969875865D+01
	  w(20)=0.10618816874131427052674325750558637D+01
	  w(21)=0.10881609252151220879056836565723643D+01
	  w(22)=0.11286083536458515474976138648344204D+01
	endif
	if(n.eq. 23) then
	  x(1)=0.41972147681820962117349438358360203D-01
	  x(2)=0.21672584044775683453920578358520288D+00
	  x(3)=0.51594567024719162340318234782898799D+00
	  x(4)=0.92153380268569184963212137517984957D+00
	  x(5)=0.14166238549372345278953844172065146D+01
	  x(6)=0.19872569198613749819060103618539870D+01
	  x(7)=0.26222049588938287634221069311024621D+01
	  x(8)=0.33123818279738364012132478819694370D+01
	  x(9)=0.40503474579712408498721742944326252D+01
	  x(10)=0.48299708035495152458338452561292797D+01
	  x(11)=0.56462005126847056857286588638089597D+01
	  x(12)=0.64948944390679841376368131022900343D+01
	  x(13)=0.73726825911123006207503749465104192D+01
	  x(14)=0.82768541685341041613810375565662980D+01
	  x(15)=0.92052664234629180128877123934216797D+01
	  x(16)=0.10156275473759043492805176356341690D+02
	  x(17)=0.11128690172219716458812399650923908D+02
	  x(18)=0.12121751229071891842181685206014663D+02
	  x(19)=0.13135139174255273886160466645378619D+02
	  x(20)=0.14169014559581503576168870495166630D+02
	  x(21)=0.15224088692195556049568949674721807D+02
	  x(22)=0.16301795992929449141684017376974225D+02
	  x(23)=0.17406564924714970743480080273002386D+02
	  w(1)=0.10717641518184028992521916734403931D+00
	  w(2)=0.23990532697980757048661359931429615D+00
	  w(3)=0.35542521636035206977766165437060408D+00
	  w(4)=0.45291831233090568131771647131245118D+00
	  w(5)=0.53494156757723620998490332567598671D+00
	  w(6)=0.60446400516340048625352210365235806D+00
	  w(7)=0.66392601599188527661254966005799361D+00
	  w(8)=0.71519317335927368173287277386407368D+00
	  w(9)=0.75971995941075953773236051347339526D+00
	  w(10)=0.79868765250436368763331529407878406D+00
	  w(11)=0.83308373888879982960276038284064271D+00
	  w(12)=0.86374448602560738308397958462592214D+00
	  w(13)=0.89138192927548665078774092762614600D+00
	  w(14)=0.91660624149239344760076164675410837D+00
	  w(15)=0.93994758288643587285093872196739540D+00
	  w(16)=0.96187891563469807287845014798222110D+00
	  w(17)=0.98284121163181492608629241658491082D+00
	  w(18)=0.10032739246903483731898631958756596D+01
	  w(19)=0.10236564911929821874991830554790795D+01
	  w(20)=0.10445723165309745983364564381190576D+01
	  w(21)=0.10668357475893164298952342505799606D+01
	  w(22)=0.10921225216001932967202492363867350D+01
	  w(23)=0.11311366331613983415849133962183259D+01
	endif
	if(n.eq. 24) then
	  x(1)=0.40161660870322704597334961817978183D-01
	  x(2)=0.20770303300921480471608049356291303D+00
	  x(3)=0.49556233372164454342723161062167492D+00
	  x(4)=0.88721719544850086158760404941858724D+00
	  x(5)=0.13669313729767131881942532345419750D+01
	  x(6)=0.19214618684060367215238329663407058D+01
	  x(7)=0.25400363129156677821640641923295334D+01
	  x(8)=0.32138670492582850357621282455511391D+01
	  x(9)=0.39357043755757925573846023326041177D+01
	  x(10)=0.46995238977378539502183230069931597D+01
	  x(11)=0.55003147577352873298650592914782465D+01
	  x(12)=0.63339248942996526636761700501665473D+01
	  x(13)=0.71969374920866178158007642196025699D+01
	  x(14)=0.80865680629338196894195789354853332D+01
	  x(15)=0.90005792910157271791149469208903611D+01
	  x(16)=0.99372135842936462069019398768432438D+01
	  x(17)=0.10895144187912299216236533538904041D+02
	  x(18)=0.11873446478870155473828162939753383D+02
	  x(19)=0.12871592488448108326792862499132752D+02
	  x(20)=0.13889473955101772517650715599302202D+02
	  x(21)=0.14927461287658761790453354478813708D+02
	  x(22)=0.15986507931569594376242093858309090D+02
	  x(23)=0.17068406196932564000690035754814744D+02
	  x(24)=0.18178282867325837202088223421014845D+02
	  w(1)=0.10259513445157118771700055503970361D+00
	  w(2)=0.23033172049800118563744888433575397D+00
	  w(3)=0.34254832958013686772247297085414175D+00
	  w(4)=0.43811525059456785102796061437402386D+00
	  w(5)=0.51910931505333612889785399602260441D+00
	  w(6)=0.58816539155050029386018195509677753D+00
	  w(7)=0.64752645162427036051866480192984454D+00
	  w(8)=0.69893183105287592926657680436619557D+00
	  w(9)=0.74374217828666766560985479372902773D+00
	  w(10)=0.78306425368614274162837318726815283D+00
	  w(11)=0.81782751820985366553173889769823290D+00
	  w(12)=0.84882505534139573999397043735370971D+00
	  w(13)=0.87673762348836203894819618653855287D+00
	  w(14)=0.90215211847778098430694626586046070D+00
	  w(15)=0.92557913449000162042779038529261015D+00
	  w(16)=0.94747115349876975809451096210977994D+00
	  w(17)=0.96824218146576768173616756030241959D+00
	  w(18)=0.98829035217221483655691827152622864D+00
	  w(19)=0.10080268468618629462696389964548871D+01
	  w(20)=0.10279180200015789470313620768138207D+01
	  w(21)=0.10485554182563696645758000158821233D+01
	  w(22)=0.10708053477833299194088567674043588D+01
	  w(23)=0.10965139925655209562904701670049690D+01
	  w(24)=0.11368396001481675394728654282516800D+01
	endif
	if(n.eq. 25) then
	  x(1)=0.38544654960740845628386352927918779D-01
	  x(2)=0.19964278415248942755155781014764216D+00
	  x(3)=0.47736849737003589355666122173715848D+00
	  x(4)=0.85665497240449139937368272512685508D+00
	  x(5)=0.13228177791134794194505275299889036D+01
	  x(6)=0.18632649417618365461635221436154097D+01
	  x(7)=0.24676211718997795863117517001228407D+01
	  x(8)=0.31273523012782726837599511782173067D+01
	  x(9)=0.38353727466903011489307573356200010D+01
	  x(10)=0.45857556281779983109458953549619764D+01
	  x(11)=0.53735328549237166839702695142477751D+01
	  x(12)=0.61945499685463563110943141509778798D+01
	  x(13)=0.70453513302110204463701847998891026D+01
	  x(14)=0.79230842038679858418959156551863998D+01
	  x(15)=0.88254179188323949034611359820701182D+01
	  x(16)=0.97504775019799900803718628594651818D+01
	  x(17)=0.10696792232004669997991186392027885D+02
	  x(18)=0.11663260073103119296433760609943420D+02
	  x(19)=0.12649129713037108757589521701447666D+02
	  x(20)=0.13654003334795369539733655983582139D+02
	  x(21)=0.14677865231928679534689763386268169D+02
	  x(22)=0.15721142576564773207792313769459724D+02
	  x(23)=0.16784804136239745986358684604056180D+02
	  x(24)=0.17870594977045829665485143777914345D+02
	  x(25)=0.18983553683850214355288699152879417D+02
	  w(1)=0.98502860786019150185488513216114370D-01
	  w(2)=0.22178069453389132048570786537311506D+00
	  w(3)=0.33108283188094378424182195885805413D+00
	  w(4)=0.42501813294363627315775033821410034D+00
	  w(5)=0.50520938140300153396111682013724931D+00
	  w(6)=0.57396300243312403654982745138113387D+00
	  w(7)=0.63333383047793623266130680349306203D+00
	  w(8)=0.68495226078388837809285405455739237D+00
	  w(9)=0.73010405660473898237228240759577602D+00
	  w(10)=0.76983605703678692488978185792802833D+00
	  w(11)=0.80502788484662313184259119225316681D+00
	  w(12)=0.83643201588340776453378566657193005D+00
	  w(13)=0.86469665126457539994220269363722764D+00
	  w(14)=0.89038188933874184982641963870264590D+00
	  w(15)=0.91397423285270529280666096383356489D+00
	  w(16)=0.93590124596355950270520906997262500D+00
	  w(17)=0.95654701404883446524962664625491016D+00
	  w(18)=0.97626912139202859375330945113091730D+00
	  w(19)=0.99541868057412585812926408834755421D+00
	  w(20)=0.10143665935320542725861514554708265D+01
	  w(21)=0.10335422972866488144916274904971942D+01
	  w(22)=0.10534977717277278497221004727180116D+01
	  w(23)=0.10750429943935644327979161971597932D+01
	  w(24)=0.10999239540083021182681477512232959D+01
	  w(25)=0.11393286143253755771098667537444271D+01
	endif
	if(n.eq. 26) then
	  x(1)=0.37057019538159492733075950354759698D-01
	  x(2)=0.19219683859954905802425173533265479D+00
	  x(3)=0.46045971214896902479196683088957798D+00
	  x(4)=0.82805130101422141919442765356507152D+00
	  x(5)=0.12812122994478714055333057331154123D+01
	  x(6)=0.18079201927629728885449367226101458D+01
	  x(7)=0.23981472807433310556746164365904406D+01
	  x(8)=0.30435901230658175542487242637434974D+01
	  x(9)=0.37373274292409601748943259735824540D+01
	  x(10)=0.44735476894021157789893550216220319D+01
	  x(11)=0.52473551816946653048034931998699903D+01
	  x(12)=0.60546294862094365285543062782380730D+01
	  x(13)=0.68919164879597207828965110820718110D+01
	  x(14)=0.77563386070883817424714834487531334D+01
	  x(15)=0.86455191519599416949404258048161864D+01
	  x(16)=0.95575192961392403390163963194936514D+01
	  x(17)=0.10490787606167046419614052865654230D+02
	  x(18)=0.11444122623412685157973101013340056D+02
	  x(19)=0.12416649553950447781858201778959483D+02
	  x(20)=0.13407813117883241815775363647844642D+02
	  x(21)=0.14417390388944724577413580846041441D+02
	  x(22)=0.15445530168678835991613595979288220D+02
	  x(23)=0.16492828612411702948747915797866881D+02
	  x(24)=0.17560456489260992896106472471728921D+02
	  x(25)=0.18650464841062735388277360470965505D+02
	  x(26)=0.19768476866194163221734925173223019D+02
	  w(1)=0.94733963379001084881281258276430890D-01
	  w(2)=0.21384206006426126345232319181377534D+00
	  w(3)=0.32031528543989479240750029021000955D+00
	  w(4)=0.41254929390710237679940064481343143D+00
	  w(5)=0.49176691815620643399853406663169153D+00
	  w(6)=0.55998309037173654711949666307191364D+00
	  w(7)=0.61909314036708051887103465560358018D+00
	  w(8)=0.67064351982741110536068163128220476D+00
	  w(9)=0.71586567032066061866402151281363331D+00
	  w(10)=0.75576118553096371943667008963529952D+00
	  w(11)=0.79116885492294941606417069124290720D+00
	  w(12)=0.82280556212477184985232270264532417D+00
	  w(13)=0.85129012269432990933637483976781368D+00
	  w(14)=0.87715909928110047388827297254465520D+00
	  w(15)=0.90087981520398463519683218692080118D+00
	  w(16)=0.92286282936149244271462066535605118D+00
	  w(17)=0.94347471535979188228537850591237657D+00
	  w(18)=0.96305166489155957787460238250787370D+00
	  w(19)=0.98191478773736651231018868202227168D+00
	  w(20)=0.10003889128129097585429008177015930D+01
	  w(21)=0.10188284918868633077693175437161699D+01
	  w(22)=0.10376578150755357832224490266526118D+01
	  w(23)=0.10574411346568342029428322348394431D+01
	  w(24)=0.10790382469712218593826946744229645D+01
	  w(25)=0.11043433786820813224949233699589968D+01
	  w(26)=0.11448816650689648266592257641605102D+01
	endif
	if(n.eq. 27) then
	  x(1)=0.35834993247954083361861421508365311D-01
	  x(2)=0.18605352495124757861155728733137948D+00
	  x(3)=0.44642194330546725034025712375296280D+00
	  x(4)=0.80414605853435716653621057048439980D+00
	  x(5)=0.12462239924414626468518463298096322D+01
	  x(6)=0.17611257461239471222569363817456178D+01
	  x(7)=0.23391453785205094106913747964426875D+01
	  x(8)=0.29722080942457780317056403873721138D+01
	  x(9)=0.36535650947301854252202701900387183D+01
	  x(10)=0.43775326951870896508012265258003026D+01
	  x(11)=0.51393010711012010460763121955096722D+01
	  x(12)=0.59347963207816585295972799940500408D+01
	  x(13)=0.67605758319435445002909546019509435D+01
	  x(14)=0.76137438021356000916739503736607730D+01
	  x(15)=0.84918804165657331139982488821260631D+01
	  x(16)=0.93929821745587691594892021385021508D+01
	  x(17)=0.10315412711899460518338855763431638D+02
	  x(18)=0.11257864229560141211550217121839523D+02
	  x(19)=0.12219330080395609527954547957051545D+02
	  x(20)=0.13199089643751765521528795943595469D+02
	  x(21)=0.14196707628360135444722800457384437D+02
	  x(22)=0.15212051628819105886236684455070645D+02
	  x(23)=0.16245334166526408381514556822367013D+02
	  x(24)=0.17297187705846923222452460322529078D+02
	  x(25)=0.18368783996172027173088281415402889D+02
	  x(26)=0.19462114300213894324542707181535661D+02
	  x(27)=0.20582673852051787122263704077340662D+02
	  w(1)=0.91634452557743162337544617912499234D-01
	  w(2)=0.20725798653776011937210910218709614D+00
	  w(3)=0.31128301536359448409285732850548811D+00
	  w(4)=0.40197254484414085551335915624804329D+00
	  w(5)=0.48025776040426659818649568478576839D+00
	  w(6)=0.54793227641803410321585943165700883D+00
	  w(7)=0.60676658418922713167376059573143721D+00
	  w(8)=0.65823842224192508076896501734154299D+00
	  w(9)=0.70353143374483517735740178977721371D+00
	  w(10)=0.74360437862720618618084245099453256D+00
	  w(11)=0.77925510412986809249957786960294470D+00
	  w(12)=0.81116302338313539532776985652162693D+00
	  w(13)=0.83991454783642205406124503497267142D+00
	  w(14)=0.86601910131774351153666202662861906D+00
	  w(15)=0.88992107452719526605733335600234568D+00
	  w(16)=0.91201046747319902685546821885509416D+00
	  w(17)=0.93263337746385499560375365035724826D+00
	  w(18)=0.95210285562789409574691035231808200D+00
	  w(19)=0.97071063776060639849418976154993288D+00
	  w(20)=0.98874065527360133209811010601697490D+00
	  w(21)=0.10064860790234053666836189222522080D+01
	  w(22)=0.10242732841383590525197178067173809D+01
	  w(23)=0.10424992963077313046227345694205724D+01
	  w(24)=0.10616963916554709168593717549811117D+01
	  w(25)=0.10826735650889751649827985602314584D+01
	  w(26)=0.11072236156410122376314575376454741D+01
	  w(27)=0.11468295974452529240039666547090746D+01
	endif
	if(n.eq. 28) then
	  x(1)=0.34649421149527286079017329711859929D-01
	  x(2)=0.18005833416581945582990442744630855D+00
	  x(3)=0.43260408861564048077497091071563773D+00
	  x(4)=0.78038984122852594893515743024181575D+00
	  x(5)=0.12111271013211080127547347728977911D+01
	  x(6)=0.17137747630321893410609845886938274D+01
	  x(7)=0.22789600476402123163666146865580231D+01
	  x(8)=0.28988505738640957254403929255204275D+01
	  x(9)=0.35668835527411730801361500198254362D+01
	  x(10)=0.42775174285095580728466302389279008D+01
	  x(11)=0.50260444316599945935308824118692428D+01
	  x(12)=0.58084554651427726668089235317893326D+01
	  x(13)=0.66213397270308567144070366339292377D+01
	  x(14)=0.74618055548714492530848474416416138D+01
	  x(15)=0.83274148871323045995040956768207252D+01
	  x(16)=0.92161279178902919539950744365341961D+01
	  x(17)=0.10126256754412748861682302958797663D+02
	  x(18)=0.11056427893157508179911019396968186D+02
	  x(19)=0.12005553805399065581127615587320179D+02
	  x(20)=0.12972814335375019823004549834877253D+02
	  x(21)=0.13957649262705848158816479553934187D+02
	  x(22)=0.14959764521043945606493252853397280D+02
	  x(23)=0.15979156486798860470344152417965233D+02
	  x(24)=0.17016161724139383437659489572979510D+02
	  x(25)=0.18071543406684014598795329220592976D+02
	  x(26)=0.19146633395269574862140871118754148D+02
	  x(27)=0.20243668518906520858990916167385876D+02
	  x(28)=0.21368602329962925523432204499840736D+02
	  w(1)=0.88622857693915105259385711633512983D-01
	  w(2)=0.20078714572495373613669755741284462D+00
	  w(3)=0.30226336710228968174263286528002936D+00
	  w(4)=0.39123451482265819700856468443816993D+00
	  w(5)=0.46839451945520893527685757362633012D+00
	  w(6)=0.53534173858832079240244183893082663D+00
	  w(7)=0.59372828885713024860137920768465847D+00
	  w(8)=0.64496316544847820395602866483386606D+00
	  w(9)=0.69018228980436679975696279143448919D+00
	  w(10)=0.73030303949562380694260355085134506D+00
	  w(11)=0.76608442505739382522023106503183953D+00
	  w(12)=0.79817047815891095652318654174450785D+00
	  w(13)=0.82711718964233071904601501955767162D+00
	  w(14)=0.85340904918321747185672165869618766D+00
	  w(15)=0.87747043602075713852883609433774836D+00
	  w(16)=0.89967497620477510089642692037159577D+00
	  w(17)=0.92035435534202525431624053453560919D+00
	  w(18)=0.93980723716479275253732339479029179D+00
	  w(19)=0.95830868989922746159493271989049390D+00
	  w(20)=0.97612065905778322427721604981343262D+00
	  w(21)=0.99350448475457564523338760409387760D+00
	  w(22)=0.10107373633334020279761489291558973D+01
	  w(23)=0.10281364360012410497091650540824048D+01
	  w(24)=0.10460977891287441821077663917094469D+01
	  w(25)=0.10651662838357145446366303076501936D+01
	  w(26)=0.10861931437449714454857030432322063D+01
	  w(27)=0.11110870028665837949688466324005276D+01
	  w(28)=0.11515563125174608227752059974591248D+01
	endif
	if(n.eq. 29) then
	  x(1)=0.33553959839036001655543373090040404D-01
	  x(2)=0.17450396418281238775094266202358995D+00
	  x(3)=0.41975868958520995422034616240125615D+00
	  x(4)=0.75824229731656256348060196614824235D+00
	  x(5)=0.11783521199560875736267462343676016D+01
	  x(6)=0.16695400125587793738191066950093955D+01
	  x(7)=0.22227783087991492827484307781560346D+01
	  x(8)=0.28304848302713394581076045142253861D+01
	  x(9)=0.34862844007039122118385421345010400D+01
	  x(10)=0.41847751598255813476612274826038629D+01
	  x(11)=0.49213473609559184041017942945472896D+01
	  x(12)=0.56920524484299326317682243825402111D+01
	  x(13)=0.64935072997818688023130562214646488D+01
	  x(14)=0.73228202200078840178321115672588348D+01
	  x(15)=0.81775303561968240018131837132386863D+01
	  x(16)=0.90555563983894327151347170001827180D+01
	  x(17)=0.99551529304501737982491249567829072D+01
	  x(18)=0.10874874011268980567024300398770720D+02
	  x(19)=0.11813544104522151201308588497340679D+02
	  x(20)=0.12770236786942932738497802347410470D+02
	  x(21)=0.13744262037156889988409602665342391D+02
	  x(22)=0.14735163540289747174938383977860212D+02
	  x(23)=0.15742728540739239662116233375854790D+02
	  x(24)=0.16767014578429062510167568689212203D+02
	  x(25)=0.17808400110182436293371210922487080D+02
	  x(26)=0.18867668960609286443741439143195748D+02
	  x(27)=0.19946143700718788949188819970004261D+02
	  x(28)=0.21045995786123555149060848634690046D+02
	  x(29)=0.22173033865751079929395928047597408D+02
	  w(1)=0.85838158802988701046032815611397382D-01
	  w(2)=0.19477385901158245085085241043998394D+00
	  w(3)=0.29383621766506690109466148896899540D+00
	  w(4)=0.38117224372587515945909331094298977D+00
	  w(5)=0.45728310074921790784330255519307684D+00
	  w(6)=0.52359020598427119530526852031471208D+00
	  w(7)=0.58162702319516323967718562926165760D+00
	  w(8)=0.63272747067482759675272063759621233D+00
	  w(9)=0.67797436135679445268209519781521522D+00
	  w(10)=0.71824100693600445577402524577337317D+00
	  w(11)=0.75424686608660118647406989111914299D+00
	  w(12)=0.78660066973655529132258834579261020D+00
	  w(13)=0.81582808550718999285322752257343382D+00
	  w(14)=0.84238846687921586919856054009869695D+00
	  w(15)=0.86668565323448476522827377266366966D+00
	  w(16)=0.88907615047281940423573587395367213D+00
	  w(17)=0.90987644549942581218004988841130398D+00
	  w(18)=0.92937024251990218992602876824093983D+00
	  w(19)=0.94781599269424532394623383879661560D+00
	  w(20)=0.96545503882224525948885229809093289D+00
	  w(21)=0.98252089587658952840598658440285362D+00
	  w(22)=0.99925063547801429297834374665399082D+00
	  w(23)=0.10159001898755433135335124461562373D+01
	  w(24)=0.10327670289817463267212360733537935D+01
	  w(25)=0.10502268960351484050619319532415830D+01
	  w(26)=0.10687987121279187707756364034139551D+01
	  w(27)=0.10892886965844568170069806001265533D+01
	  w(28)=0.11135100571950677927901551811373793D+01
	  w(29)=0.11530882719782502565664117355481721D+01
	endif
	if(n.eq. 30) then
	  x(1)=0.32395423845231328952021243594572297D-01
	  x(2)=0.16861844033714026847015077237301739D+00
	  x(3)=0.40611377169029516664977563777938485D+00
	  x(4)=0.73466473057596815010583668481558561D+00
	  x(5)=0.11434056199839885969993247272213921D+01
	  x(6)=0.16223240841225237129918923528748564D+01
	  x(7)=0.21627613886742245163929965201532468D+01
	  x(8)=0.27573919900368202284823837544536218D+01
	  x(9)=0.34000247011207842540159163036150858D+01
	  x(10)=0.40853910479355226570419290510471910D+01
	  x(11)=0.48089751549709580658031882194336504D+01
	  x(12)=0.55668891598344441007384375552646816D+01
	  x(13)=0.63557824365416699663455801783129573D+01
	  x(14)=0.71727723299071382001557140029035509D+01
	  x(15)=0.80153880354211217706961178919300437D+01
	  x(16)=0.88815231304950224000549496849998832D+01
	  x(17)=0.97693948098293734716435210430063307D+01
	  x(18)=0.10677509220347507934434361231978983D+02
	  x(19)=0.11604632899927898037617524096276611D+02
	  x(20)=0.12549770612996521634840974002145231D+02
	  x(21)=0.13512150122572977650747816369403154D+02
	  x(22)=0.14491214826551962246981020143721253D+02
	  x(23)=0.15486625876302246851423660700675100D+02
	  x(24)=0.16498276597704045087766644428484142D+02
	  x(25)=0.17526324055306258031805555219762027D+02
	  x(26)=0.18571245797007211564277895377017558D+02
	  x(27)=0.19633934281183009318283438915386796D+02
	  x(28)=0.20715851636750954867238760925829411D+02
	  x(29)=0.21819391138662258100566759821958840D+02
	  x(30)=0.22950804950088933509277921984903514D+02
	  w(1)=0.82891596110067619629369062295154436D-01
	  w(2)=0.18838810673274555829159737641020911D+00
	  w(3)=0.28485143005306101660067952252575196D+00
	  w(4)=0.37041553715895586540440831413434353D+00
	  w(5)=0.44539043894975510751166325462691020D+00
	  w(6)=0.51100452150290254227371633533039130D+00
	  w(7)=0.56865283856139159368581204034853727D+00
	  w(8)=0.61958013174010084789244956482434645D+00
	  w(9)=0.66481004321965575520891889027552679D+00
	  w(10)=0.70517204769960306620646406372543424D+00
	  w(11)=0.74134967169016696697525503623182885D+00
	  w(12)=0.77392103530415079504223285766784102D+00
	  w(13)=0.80338600122756931298795279872138053D+00
	  w(14)=0.83018277269650053717953142040641978D+00
	  w(15)=0.85469824839953001394121656630886719D+00
	  w(16)=0.87727539085565742205119477148400620D+00
	  w(17)=0.89821948245755822703983994870213792D+00
	  w(18)=0.91780416582368185007823058185749687D+00
	  w(19)=0.93627766216629293438700187834911048D+00
	  w(20)=0.95386940504388062489482535966089927D+00
	  w(21)=0.97079739700556766379690998292062432D+00
	  w(22)=0.98727684670885629270031813575769775D+00
	  w(23)=0.10035311243345927945114226531586610D+01
	  w(24)=0.10198069790571226977959895521053113D+01
	  w(25)=0.10363977445722241110814820785890333D+01
	  w(26)=0.10536819126632297471957144807674922D+01
	  w(27)=0.10721934390392973135419651953270659D+01
	  w(28)=0.10927831816201418746459239628165960D+01
	  w(29)=0.11173737370677914171324118797201663D+01
	  w(30)=0.11578618493114158916057476744754240D+01
	endif
	if(n.eq. 31) then
	  x(1)=0.31354278310343070002019771891355049D-01
	  x(2)=0.16332835712339532907577677178778686D+00
	  x(3)=0.39385199392482811431293043824553024D+00
	  x(4)=0.71349775214722199745409625393222086D+00
	  x(5)=0.11120868650386663833984357552253641D+01
	  x(6)=0.15801071341254320223868035100167617D+01
	  x(7)=0.21092378289133741908756292104953900D+01
	  x(8)=0.26923802662837172938736784999491647D+01
	  x(9)=0.33234998137128841655396627174923196D+01
	  x(10)=0.39974385925377480255349382787244394D+01
	  x(11)=0.47097577456400570028449692472349852D+01
	  x(12)=0.54566190267492231669166358187794685D+01
	  x(13)=0.62346978423822019266253846581093967D+01
	  x(14)=0.70411172062729967535688047064468265D+01
	  x(15)=0.78733946685933355880138151405844837D+01
	  x(16)=0.87293975824146112785228979191742837D+01
	  x(17)=0.96073044951270674118859460577368736D+01
	  x(18)=0.10505571845959279286830678756814450D+02
	  x(19)=0.11422905824103905203514841559808701D+02
	  x(20)=0.12358239542725163673253518936689943D+02
	  x(21)=0.13310715865946766811589441203977913D+02
	  x(22)=0.14279676452924514151732182654086500D+02
	  x(23)=0.15264657976271370998233578575309366D+02
	  x(24)=0.16265397160565473910764922038652003D+02
	  x(25)=0.17281847462814127425190235953778028D+02
	  x(26)=0.18314212155925179814630610053427517D+02
	  x(27)=0.19363001659818237953913921955972910D+02
	  x(28)=0.20429126767095223016212912625633180D+02
	  x(29)=0.21514047061292600915294315200299025D+02
	  x(30)=0.22620126001584271335786979761905968D+02
	  x(31)=0.23753604827909729380053249769844115D+02
	  w(1)=0.80243398875130558800528035590104992D-01
	  w(2)=0.18264806785017626000744428438338218D+00
	  w(3)=0.27678431664690572588938266562763602D+00
	  w(4)=0.36078787836066467109574773530766834D+00
	  w(5)=0.43479728577654619403602964666788466D+00
	  w(6)=0.49985669209612632046813018860120792D+00
	  w(7)=0.55722457082888166013390218722634017D+00
	  w(8)=0.60805777233524505032136175941559486D+00
	  w(9)=0.65332431619085751250963767233770341D+00
	  w(10)=0.69381538132515707673064753180369735D+00
	  w(11)=0.73018463967704683170723001239821315D+00
	  w(12)=0.76298490548565767177535690279910341D+00
	  w(13)=0.79269410213246049856650188303319737D+00
	  w(14)=0.81973162509324370983421204073238187D+00
	  w(15)=0.84446850104307125128855204820865765D+00
	  w(16)=0.86723430801843581772203606305993162D+00
	  w(17)=0.88832273165677944071205729414941743D+00
	  w(18)=0.90799672835074207721817174387979321D+00
	  w(19)=0.92649373188643324716906590765574947D+00
	  w(20)=0.94403111330718736393663448325241916D+00
	  w(21)=0.96081208209329649339025536392000504D+00
	  w(22)=0.97703233520897692088169605995062739D+00
	  w(23)=0.99288801493529488251255088471225463D+00
	  w(24)=0.10085859969498260912956766333081760D+01
	  w(25)=0.10243583939923772696545256621902809D+01
	  w(26)=0.10404848104794099938885665324050933D+01
	  w(27)=0.10573296383546104681983024420333095D+01
	  w(28)=0.10754089559129176301155439432477579D+01
	  w(29)=0.10955423446690406574788312354939990D+01
	  w(30)=0.11196138285621031194239094475051388D+01
	  w(31)=0.11596321256631261231717644477612339D+01
	endif
	if(n.eq. 32) then
	  x(1)=0.30349937742381039029737621603999287D-01
	  x(2)=0.15821530524930957350981941544887377D+00
	  x(3)=0.38196519764129310825495622339076363D+00
	  x(4)=0.69290532029415574655217824329156429D+00
	  x(5)=0.10815000835196373341062781037180685D+01
	  x(6)=0.15387055890220047160710237221792340D+01
	  x(7)=0.20565151130412311353268250968540087D+01
	  x(8)=0.26280393793364540400148143817204982D+01
	  x(9)=0.32473899625462565587952212808886543D+01
	  x(10)=0.39095168843281236092934705084189773D+01
	  x(11)=0.46100629001653583571851413580588996D+01
	  x(12)=0.53452500349782114952290612563956529D+01
	  x(13)=0.61117948249824189943524288537446409D+01
	  x(14)=0.69068439706878734796191565692424774D+01
	  x(15)=0.77279231924730771297049614076968282D+01
	  x(16)=0.85728946088327564467590491403825581D+01
	  x(17)=0.94399201347038257381427683867514133D+01
	  x(18)=0.10327429796843327736155515594873577D+02
	  x(19)=0.11234094613649570959523771307431161D+02
	  x(20)=0.12158804045203877919334445323329419D+02
	  x(21)=0.13100648205518305999817130214069039D+02
	  x(22)=0.14058905194765479507168493000790477D+02
	  x(23)=0.15033034160318033300995921308640391D+02
	  x(24)=0.16022675127681903717302702716551721D+02
	  x(25)=0.17027657384458382949787846882827580D+02
	  x(26)=0.18048019461323221435122832190245390D+02
	  x(27)=0.19084045896883136350652421242557466D+02
	  x(28)=0.20136329387464179774269723566249013D+02
	  x(29)=0.21205872248485881925716967089101672D+02
	  x(30)=0.22294251962151093948705238290131092D+02
	  x(31)=0.23404022785632438541369992890395224D+02
	  x(32)=0.24541777752285696578837814740836620D+02
	  w(1)=0.77687554219704149649494695495377528D-01
	  w(2)=0.17708703705627543145695312887255568D+00
	  w(3)=0.26892445951139809512397960133966990D+00
	  w(4)=0.35134485491478623098160483095853124D+00
	  w(5)=0.42432912048747656319491738940996584D+00
	  w(6)=0.48874675549663260110122564583434723D+00
	  w(7)=0.54572636490278558785149698451277800D+00
	  w(8)=0.59634360129512142822960640842211433D+00
	  w(9)=0.64151835689239933291361239753314294D+00
	  w(10)=0.68201115660547106944733286582049914D+00
	  w(11)=0.71845260343410932701857518623000942D+00
	  w(12)=0.75137547715096031364367945570847951D+00
	  w(13)=0.78123939449719914662750852585304528D+00
	  w(14)=0.80844722185904593914074212079867721D+00
	  w(15)=0.83335551610918434661812170816119760D+00
	  w(16)=0.85628149724405755449652133393101394D+00
	  w(17)=0.87750834620079276593429540298529901D+00
	  w(18)=0.89728986194176141300005156153929420D+00
	  w(19)=0.91585499021399385810582316480576992D+00
	  w(20)=0.93341246626577922818057686527026817D+00
	  w(21)=0.95015572913875601113176116996328346D+00
	  w(22)=0.96626830574334088286292399061494507D+00
	  w(23)=0.98193000412019415534814470447599888D+00
	  w(24)=0.99732451963323776045200474982266314D+00
	  w(25)=0.10126495345374315615316618277574889D+01
	  w(26)=0.10281312953198651349850933911511675D+01
	  w(27)=0.10440474832830413109974188046180643D+01
	  w(28)=0.10607663285146948783932430160348304D+01
	  w(29)=0.10788165311873763041461415923549794D+01
	  w(30)=0.10990633062876840053689875276177190D+01
	  w(31)=0.11234609942681621497229116357630119D+01
	  w(32)=0.11642028453735175919092625917983241D+01
	endif
	if(n.eq. 33) then
	  x(1)=0.29516841548121147686956078359798994D-01
	  x(2)=0.15396436335903893044552148694492644D+00
	  x(3)=0.37204808347099149745673685174551792D+00
	  x(4)=0.67565524712779223559522279174416326D+00
	  x(5)=0.10557709540672579606734871049411595D+01
	  x(6)=0.15037398176183069953282256392412819D+01
	  x(7)=0.20118235591875479606471799343125895D+01
	  x(8)=0.25733186070240368081840642844326794D+01
	  x(9)=0.31824715429064855953811274957843125D+01
	  x(10)=0.38343359509698453635451187437865883D+01
	  x(11)=0.45246356468237429737655475037172437D+01
	  x(12)=0.52496544518904109111190336989238858D+01
	  x(13)=0.60061525587093376188363436085637659D+01
	  x(14)=0.67913031119156554993310237477999181D+01
	  x(15)=0.76026424249740527372409815143328160D+01
	  x(16)=0.84380290860979751954573657712899148D+01
	  x(17)=0.92956091443113493255623325239866972D+01
	  x(18)=0.10173785967015565034898827434517443D+02
	  x(19)=0.11071194190674610524638410424813628D+02
	  x(20)=0.11986677616627256171000226459000260D+02
	  x(21)=0.12919271119546655413046210014726967D+02
	  x(22)=0.13868186777048956415114844276104122D+02
	  x(23)=0.14832804591874367616810559411533177D+02
	  x(24)=0.15812668437981050928442527947481722D+02
	  x(25)=0.16807488298842780949371444876305759D+02
	  x(26)=0.17817150602982401608187501551583409D+02
	  x(27)=0.18841739696439972817643138114362955D+02
	  x(28)=0.19881575591888250187366793397814035D+02
	  x(29)=0.20937275961922320988151113851927221D+02
	  x(30)=0.22009856525751274602953344583511353D+02
	  x(31)=0.23100893760260753140300948871299624D+02
	  x(32)=0.24212906935700107879938514088280499D+02
	  x(33)=0.25352458767932045446968913893215358D+02
	  w(1)=0.75566274379686090845176238417479908D-01
	  w(2)=0.17245090185343961430852743887953693D+00
	  w(3)=0.26232862929412148478292010622681119D+00
	  w(4)=0.34336258354556015159886328547145240D+00
	  w(5)=0.41541767735997026100847051566233858D+00
	  w(6)=0.47922896601846382003131452620436903D+00
	  w(7)=0.53582299405173028361559772747568786D+00
	  w(8)=0.58621072438652011005189024217543192D+00
	  w(9)=0.63127306256395476768972230274812318D+00
	  w(10)=0.67174488823670497339435314643196762D+00
	  w(11)=0.70823651688728472830547389094135724D+00
	  w(12)=0.74126216601389005589339831203687936D+00
	  w(13)=0.77126372355325145679927345554460771D+00
	  w(14)=0.79862744362095428041214972836314701D+00
	  w(15)=0.82369485776729867154699604725465178D+00
	  w(16)=0.84676997408885934515154758628341369D+00
	  w(17)=0.86812447146877491643834900969523005D+00
	  w(18)=0.88800198359922455360049298178637400D+00
	  w(19)=0.90662207006592154812096850946545601D+00
	  w(20)=0.92418416998312613586108454910572618D+00
	  w(21)=0.94087169792305158644296625425340608D+00
	  w(22)=0.95685642254984970556108692107954994D+00
	  w(23)=0.97230333230250642895242663144017570D+00
	  w(24)=0.98737633584223760685461002140073106D+00
	  w(25)=0.10022453993514741643622301126015373D+01
	  w(26)=0.10170961830507487810848488152259961D+01
	  w(27)=0.10321441128244432938743102567968890D+01
	  w(28)=0.10476567228822977195790144833154045D+01
	  w(29)=0.10639899573075621841411475543281995D+01
	  w(30)=0.10816619934110704992491491793771274D+01
	  w(31)=0.11014955209885703979466597957070917D+01
	  w(32)=0.11254225677442042385933973491773941D+01
	  w(33)=0.11656987034000894087881761151948012D+01
	endif
	if(n.eq. 34) then
	  x(1)=0.28692651028194224305245541017939104D-01
	  x(2)=0.14974476358725391511583779902139213D+00
	  x(3)=0.36215242246831369499204811290837824D+00
	  x(4)=0.65833467639295095708007465873379260D+00
	  x(5)=0.10297667251430566714986980514368042D+01
	  x(6)=0.14681718575272575399992547318106517D+01
	  x(7)=0.19660836342145164845618410254246555D+01
	  x(8)=0.25169935589275236331729956873459741D+01
	  x(9)=0.31152944132230948603989872935926542D+01
	  x(10)=0.37561547919856645982861209631664678D+01
	  x(11)=0.44353910610802671499186544679105282D+01
	  x(12)=0.51493605648020599474534719774965197D+01
	  x(13)=0.58948792661202604747927580319810659D+01
	  x(14)=0.66691592755782664880825905129313469D+01
	  x(15)=0.74697603913674583964166231453418732D+01
	  x(16)=0.82945509332970335236723258276470006D+01
	  x(17)=0.91416748168513723982187002548016608D+01
	  x(18)=0.10009523178490880468416435178369284D+02
	  x(19)=0.10896709755369746730480073892977089D+02
	  x(20)=0.11802049726291521380971971666440368D+02
	  x(21)=0.12724541982253466443353318027220666D+02
	  x(22)=0.13663354942079559251055798085872084D+02
	  x(23)=0.14617816153939974910258570162113756D+02
	  x(24)=0.15587406094181691429412239813245833D+02
	  x(25)=0.16571756852068077137118962127715349D+02
	  x(26)=0.17570656852874737552383521688170731D+02
	  x(27)=0.18584063548369215368438744917511940D+02
	  x(28)=0.19612127322511650362457658047787845D+02
	  x(29)=0.20655232129740205948564835125580430D+02
	  x(30)=0.21714061832138220609067502664402127D+02
	  x(31)=0.22789708150232893757447527605108917D+02
	  x(32)=0.23883847152355446752380885300226510D+02
	  x(33)=0.24999160469212206692191102774813771D+02
	  x(34)=0.26142540164298829097333509707823396D+02
	  w(1)=0.73465920097870257565730867099773604D-01
	  w(2)=0.16783021269553297094567767544504022D+00
	  w(3)=0.25568890989410358871580797313072253D+00
	  w(4)=0.33523521193298011500161237563588656D+00
	  w(5)=0.40624063952743300065506559803907294D+00
	  w(6)=0.46932328427959341121677994124183897D+00
	  w(7)=0.52541677947477904986328667291672900D+00
	  w(8)=0.57547077013115410792920556559693068D+00
	  w(9)=0.62032788637594793446083940580138005D+00
	  w(10)=0.66069738134281064123598525839042850D+00
	  w(11)=0.69716940706570784680451424719649367D+00
	  w(12)=0.73023995388151263785658784399856813D+00
	  w(13)=0.76033359683833701136990157465334050D+00
	  w(14)=0.78782036745308980929536346593522467D+00
	  w(15)=0.81302712908770957689341685181716457D+00
	  w(16)=0.83624507863064123647234282543649897D+00
	  w(17)=0.85773495454381520897157997751492076D+00
	  w(18)=0.87773107816218653987050402065506205D+00
	  w(19)=0.89644490268434107171913183265132830D+00
	  w(20)=0.91406842840044344455918690073303878D+00
	  w(21)=0.93077767046171333564785754788317718D+00
	  w(22)=0.94673630331971581064465226518223062D+00
	  w(23)=0.96209962144106575543389681115513667D+00
	  w(24)=0.97701903437855897038843977497890592D+00
	  w(25)=0.99164746386879987216644849468139000D+00
	  w(26)=0.10061462770858133897178277038619854D+01
	  w(27)=0.10206948638656243577571558489580639D+01
	  w(28)=0.10355048992731354928764631040394306D+01
	  w(29)=0.10508431829407769431128372161765583D+01
	  w(30)=0.10670697759020835171384078421397135D+01
	  w(31)=0.10847160169351408587345986234140582D+01
	  w(32)=0.11046428228365237078634208955918439D+01
	  w(33)=0.11288634971393467498756990607944317D+01
	  w(34)=0.11697967629408208001251523455721326D+01
	endif
	if(n.eq. 35) then
	  x(1)=0.27943149453128562909043353101878893D-01
	  x(2)=0.14590069411832473633694462478160858D+00
	  x(3)=0.35311407160938601546007475917576812D+00
	  x(4)=0.64247166668018773627579776075435802D+00
	  x(5)=0.10058950842026235239501374962856062D+01
	  x(6)=0.14354652172354906980444866348989308D+01
	  x(7)=0.19239811948067189195654691502568312D+01
	  x(8)=0.24651294201415425710877116216579452D+01
	  x(9)=0.30534493384161485174388417362933978D+01
	  x(10)=0.36842234616150495440933809732086957D+01
	  x(11)=0.43533578554033427465697059233207256D+01
	  x(12)=0.50572793867785295063299599860329181D+01
	  x(13)=0.57928553206686634879929442831780761D+01
	  x(14)=0.65573322533383304744347697123885155D+01
	  x(15)=0.73482891759340853710114060959313065D+01
	  x(16)=0.81636000582367458378030278254300356D+01
	  x(17)=0.90014027461965540055643941741436720D+01
	  x(18)=0.98600722732568506501138472231104970D+01
	  x(19)=0.10738197612938781944080801622476429D+02
	  x(20)=0.11634561460544686539719805296044797D+02
	  x(21)=0.12548122934780797876896940579172224D+02
	  x(22)=0.13478003249726681289644147909712046D+02
	  x(23)=0.14423474512261989843864284921437502D+02
	  x(24)=0.15383951912487608737478694820310920D+02
	  x(25)=0.16358989743998023413951159454882145D+02
	  x(26)=0.17348281969369558908056205837056041D+02
	  x(27)=0.18351668507278066755361578543670475D+02
	  x(28)=0.19369149191722343061883293557912111D+02
	  x(29)=0.20400908642975103646222123643383384D+02
	  x(30)=0.21447357478426170018792618066072464D+02
	  x(31)=0.22509198895853941735367698129266500D+02
	  x(32)=0.23587535173489076356645455234684050D+02
	  x(33)=0.24684039873950634813581928028725088D+02
	  x(34)=0.25801405524041360450837601092644036D+02
	  x(35)=0.26946472793468675632766462513245642D+02
	  w(1)=0.71555031340840860587526606195751810D-01
	  w(2)=0.16361189172121312629393230508867418D+00
	  w(3)=0.24959951309471137270712404188088840D+00
	  w(4)=0.32775070928790778213013368258543778D+00
	  w(5)=0.39776743234904299306720076856436208D+00
	  w(6)=0.46017034363262798724036883868393488D+00
	  w(7)=0.51581021578199481858462149830302224D+00
	  w(8)=0.56557812294637710426314924916368909D+00
	  w(9)=0.61027751233960847443427155667450279D+00
	  w(10)=0.65059013286410471987863957110675983D+00
	  w(11)=0.68708446304822023975589218025561422D+00
	  w(12)=0.72023748158115508122278924929560162D+00
	  w(13)=0.75045625663456561049713400279870257D+00
	  w(14)=0.77809472328804463181484152300981805D+00
	  w(15)=0.80346529757120999093444879690650851D+00
	  w(16)=0.82684655071315404928355974334408529D+00
	  w(17)=0.84848838411988036778410560145857744D+00
	  w(18)=0.86861583591653057645487479021539912D+00
	  w(19)=0.88743224833632561310281516853137873D+00
	  w(20)=0.90512220720963532105685089845792390D+00
	  w(21)=0.92185446773691992472521405943552963D+00
	  w(22)=0.93778498662542864039437517931219190D+00
	  w(23)=0.95306016086732836534167745412560180D+00
	  w(24)=0.96782040802503321863525798107730225D+00
	  w(25)=0.98220430627346111407405260251834989D+00
	  w(26)=0.99635366284526971636381631469703279D+00
	  w(27)=0.10104201361855729146554949693381786D+01
	  w(28)=0.10245745124422205307723743317183107D+01
	  w(29)=0.10390206451950096955272329068975523D+01
	  w(30)=0.10540178410470131264986548558226787D+01
	  w(31)=0.10699188164138111467593716952251270D+01
	  w(32)=0.10872412754612261753806023989454843D+01
	  w(33)=0.11068318334853040063592288788640872D+01
	  w(34)=0.11307830968943992644426543847657740D+01
	  w(35)=0.11719202724545814131573706617928110D+01
	endif
	if(n.eq. 36) then
	  x(1)=0.27142056328972239548358302840824763D-01
	  x(2)=0.14178681907151427510349606109230081D+00
	  x(3)=0.34342381977838998263763414797722362D+00
	  x(4)=0.62543278953658920649161245819414034D+00
	  x(5)=0.98021315858273505572384465267532505D+00
	  x(6)=0.14002363417962260250249073578743264D+01
	  x(7)=0.18785945055991226659841686341678724D+01
	  x(8)=0.24091879673843461517890318646095693D+01
	  x(9)=0.29867138306149492166241543600335717D+01
	  x(10)=0.36065737319948958372606284683570266D+01
	  x(11)=0.42647654793772025172415851557161659D+01
	  x(12)=0.49577865215147900457282048591878265D+01
	  x(13)=0.56825569402033835331167210824787617D+01
	  x(14)=0.64363605312603828778605929983314127D+01
	  x(15)=0.72167996452014815389475188567303121D+01
	  x(16)=0.80217595113515507421197980875149369D+01
	  x(17)=0.88493788614634816269699513213708997D+01
	  x(18)=0.96980248504438417711526199127547443D+01
	  x(19)=0.10566271182166131481494630861561745D+02
	  x(20)=0.11452878934627758411579634412191808D+02
	  x(21)=0.12356779911330106003219952981453389D+02
	  x(22)=0.13277062520347925556052359752357006D+02
	  x(23)=0.14212960280227379783468677487690002D+02
	  x(24)=0.15163843131286117937861490645445883D+02
	  x(25)=0.16129211834035846351298459921963513D+02
	  x(26)=0.17108695899939579732063066330738366D+02
	  x(27)=0.18102055768589796258538626716472208D+02
	  x(28)=0.19109190390395511371934844646602869D+02
	  x(29)=0.20130152154913822926118882605805993D+02
	  x(30)=0.21165172450830002759403214440681040D+02
	  x(31)=0.22214703800018561707929620752111077D+02
	  x(32)=0.23279488597200625576988386455923319D+02
	  x(33)=0.24360673016746289931688806973397732D+02
	  x(34)=0.25460002587127451789683618699200451D+02
	  x(35)=0.26580225791803091084375409991480410D+02
	  x(36)=0.27728035052626292866762014455161989D+02
	  w(1)=0.69511944436498607213792411130270921D-01
	  w(2)=0.15909081406147243531457036169740604D+00
	  w(3)=0.24305240853533374711936687617708230D+00
	  w(4)=0.31968136778214911730700009684369434D+00
	  w(5)=0.38861583821018325091145584337937180D+00
	  w(6)=0.45027651490975162396068753878353164D+00
	  w(7)=0.50542366833399765546630533208372071D+00
	  w(8)=0.55488112295767544335234333630069159D+00
	  w(9)=0.59940607536473111682795433807768859D+00
	  w(10)=0.63964831820243539528547671579872258D+00
	  w(11)=0.67615269572807168430017554783262312D+00
	  w(12)=0.70937703945268570926430129475193098D+00
	  w(13)=0.73971183935845019608024131230195053D+00
	  w(14)=0.76749632153441160742346482948050834D+00
	  w(15)=0.79302992954234163835991466839914210D+00
	  w(16)=0.81658000692285714894325110435602255D+00
	  w(17)=0.83838691595559167168971725914161652D+00
	  w(18)=0.85866766530031546356127591934637167D+00
	  w(19)=0.87761878886172950409161330753704533D+00
	  w(20)=0.89541891941597706594535566182457842D+00
	  w(21)=0.91223129526075630302273111738031730D+00
	  w(22)=0.92820632453781548587556926577235572D+00
	  w(23)=0.94348428787583971111274649956612848D+00
	  w(24)=0.95819826511204764241114162359735928D+00
	  w(25)=0.97247741465258175086461278624483384D+00
	  w(26)=0.98645081746307394787720568274380639D+00
	  w(27)=0.10002522385876044808128426666371524D+01
	  w(28)=0.10140264169597386079146872361889109D+01
	  w(29)=0.10279380258279977589097597956424579D+01
	  w(30)=0.10421853786456081181199806451331824D+01
	  w(31)=0.10570235062817368021370612041209824D+01
	  w(32)=0.10728043592728884192410987452603877D+01
	  w(33)=0.10900568310850353714869243049179204D+01
	  w(34)=0.11096639729316220002175441550207324D+01
	  w(35)=0.11336107695019650432044500121264718D+01
	  w(36)=0.11734853243038561032562938635237515D+01
	endif
	if(n.eq. 37) then
	  x(1)=0.26411910931522439294205639725987567D-01
	  x(2)=0.13803646212333578469078076977893943D+00
	  x(3)=0.33458866458286440348501855623908341D+00
	  x(4)=0.60990059393603779369641415541991591D+00
	  x(5)=0.95681745451752964193303796491818503D+00
	  x(6)=0.13681782159121467223883428232511505D+01
	  x(7)=0.18373504640983484570426753634819761D+01
	  x(8)=0.23584324939925189745792977191740647D+01
	  x(9)=0.29262632101579311161287932918639854D+01
	  x(10)=0.35363476922733063645409856690093875D+01
	  x(11)=0.41847608285535615735284409311134368D+01
	  x(12)=0.48680578548268362837347922322805971D+01
	  x(13)=0.55832012213620751950315934664104134D+01
	  x(14)=0.63275038677703676626151718664914370D+01
	  x(15)=0.70985855606344072654678711842279881D+01
	  x(16)=0.78943385082265287167047063121572137D+01
	  x(17)=0.87128992260604185560168843949213624D+01
	  x(18)=0.95526246361202105106258386513218284D+01
	  x(19)=0.10412071231990294606362112972419709D+02
	  x(20)=0.11289976725642052457487807259894907D+02
	  x(21)=0.12185243941497299502429996209684759D+02
	  x(22)=0.13096926913784281865105185715947300D+02
	  x(23)=0.14024219245635096697810695332009345D+02
	  x(24)=0.14966444857664761514115525642409921D+02
	  x(25)=0.15923051329229624073491322633344680D+02
	  x(26)=0.16893606147966337971411121543496847D+02
	  x(27)=0.17877796363475876972870537429116666D+02
	  x(28)=0.18875432425143959136448756908066571D+02
	  x(29)=0.19886457443917272769340343074873090D+02
	  x(30)=0.20910963800098624432166616315953434D+02
	  x(31)=0.21949220245145550478582663345150650D+02
	  x(32)=0.23001714444389616431863032630644739D+02
	  x(33)=0.24069220006013843260461726458743215D+02
	  x(34)=0.25152902844913263891157839680090547D+02
	  x(35)=0.26254489692371379305768641643226147D+02
	  x(36)=0.27376803594474612424392034881748259D+02
	  x(37)=0.28528002198033593828085940913297236D+02
	  w(1)=0.67649676167360742184975208601827035D-01
	  w(2)=0.15496839931055050865182920460938476D+00
	  w(3)=0.23708305259617798443905201111192582D+00
	  w(4)=0.31233209049937799006713134986057412D+00
	  w(5)=0.38029950932003797614555651307455264D+00
	  w(6)=0.44131360550008746912453716504387558D+00
	  w(7)=0.49604736684752565034628446483111475D+00
	  w(8)=0.54525789333571705874703638983191922D+00
	  w(9)=0.58965605446404534273341369043919258D+00
	  w(10)=0.62986023397233292797636750037781894D+00
	  w(11)=0.66639312805720818388977022550534457D+00
	  w(12)=0.69969560255345997212828024203190580D+00
	  w(13)=0.73014400539591861605259737189044245D+00
	  w(14)=0.75806514767400290377707960942643695D+00
	  w(15)=0.78374740362428396700522625906160101D+00
	  w(16)=0.80744830649832832758505674064508639D+00
	  w(17)=0.82939962879891959701694759132806212D+00
	  w(18)=0.84981091234194716665939495214843191D+00
	  w(19)=0.86887216649901544140988107756129466D+00
	  w(20)=0.88675619025991547772491685464046896D+00
	  w(21)=0.90362077527298800827537661461974494D+00
	  w(22)=0.91961092564740731791772532233153470D+00
	  w(23)=0.93486117283187752935447178970207460D+00
	  w(24)=0.94949805240417506357886168188997544D+00
	  w(25)=0.96364283281658358326637880963971838D+00
	  w(26)=0.97741463773467263109751002048142254D+00
	  w(27)=0.99093419212349864722000347683206201D+00
	  w(28)=0.10043285584637444696198826932231896D+01
	  w(29)=0.10177374684816955507926650170702487D+01
	  w(30)=0.10313222078463066821996108046732843D+01
	  w(31)=0.10452789938153590476588306046323851D+01
	  w(32)=0.10598598197124693420789753872668371D+01
	  w(33)=0.10754097758743259571900807713973336D+01
	  w(34)=0.10924360626488454073523826082237065D+01
	  w(35)=0.11117141892961428517594413278857246D+01
	  w(36)=0.11353410551996869770619014161638916D+01
	  w(37)=0.11809043377131303032001596875488758D+01
	endif
	if(n.eq. 38) then
	  x(1)=0.25700257544476846210956821892068547D-01
	  x(2)=0.13437730671510350610731165943434462D+00
	  x(3)=0.32595440017712046243403278822370339D+00
	  x(4)=0.59469192164253537047358122435980476D+00
	  x(5)=0.93386055036037418641825524900923483D+00
	  x(6)=0.13366516749940666741736094991210848D+01
	  x(7)=0.17966966989055765946403653288143687D+01
	  x(8)=0.23082824719392154300123820576118305D+01
	  x(9)=0.28663822209047862266118045226903632D+01
	  x(10)=0.34665993632945317770577275950927287D+01
	  x(11)=0.41050836339964167365224056993611157D+01
	  x(12)=0.47784489121199582584154086362104863D+01
	  x(13)=0.54837038163375826371748189558275044D+01
	  x(14)=0.62181967441935146467812955961562693D+01
	  x(15)=0.69795731405245007294979586731642485D+01
	  x(16)=0.77657418386686707378885330399498343D+01
	  x(17)=0.85748477001384681983608970767818391D+01
	  x(18)=0.94052485653211341798396460944786668D+01
	  x(19)=0.10255495278392480429374700179323554D+02
	  x(20)=0.11124314109370851966218651796225458D+02
	  x(21)=0.12010591256336200061127783556003124D+02
	  x(22)=0.12913359319214089637739562022034079D+02
	  x(23)=0.13831785750856202099612346501089633D+02
	  x(24)=0.14765163360026418715165164030622691D+02
	  x(25)=0.15712903001068566410936000465881079D+02
	  x(26)=0.16674528663719286214472958818078041D+02
	  x(27)=0.17649675301574401231619049212895334D+02
	  x(28)=0.18638089938014424262746615568175912D+02
	  x(29)=0.19639636918084001138140592956915498D+02
	  x(30)=0.20654308688728896470365725690498948D+02
	  x(31)=0.21682244324286205738872013171203434D+02
	  x(32)=0.22723759221275248165738958050496876D+02
	  x(33)=0.23779391643480760620832370477728546D+02
	  x(34)=0.24849974804387837679087169817648828D+02
	  x(35)=0.25936753996308112846236326731741428D+02
	  x(36)=0.27041596697766728141232306370511651D+02
	  x(37)=0.28167590252521488736192623036913574D+02
	  x(38)=0.29320962649240779995807315572164953D+02
	  w(1)=0.65834105304151860349470837263652356D-01
	  w(2)=0.15094121610877292427588258760806639D+00
	  w(3)=0.23123369629337622854770017966075102D+00
	  w(4)=0.30510468405302060279637998974067159D+00
	  w(5)=0.37208937577467365853678415987815242D+00
	  w(6)=0.43242816098550262982058711713762023D+00
	  w(7)=0.48670941510737797397467829796369188D+00
	  w(8)=0.53562540884530773777072454322478734D+00
	  w(9)=0.57984278819344259314760847701109014D+00
	  w(10)=0.61995117810859579066118385526351631D+00
	  w(11)=0.65645416844019888191752443162840791D+00
	  w(12)=0.68977873388914845431685307630687021D+00
	  w(13)=0.72028984500715997896236331143882126D+00
	  w(14)=0.74830413909214332779384903915342875D+00
	  w(15)=0.77410057710162016864074985278421082D+00
	  w(16)=0.79792801073672492417898638450424187D+00
	  w(17)=0.82001035032794500345687538356287405D+00
	  w(18)=0.84055014635036973125181702926056460D+00
	  w(19)=0.85973124956521562722855378524400294D+00
	  w(20)=0.87772100427513155551650925190187991D+00
	  w(21)=0.89467224914328558060816476427135058D+00
	  w(22)=0.91072527777349121524963493357063271D+00
	  w(23)=0.92600984369793104544754669404937886D+00
	  w(24)=0.94064726902680118048039048517239280D+00
	  w(25)=0.95475272252572462861053281812928617D+00
	  w(26)=0.96843776378844259600242594387964346D+00
	  w(27)=0.98181331353225009817009549806243740D+00
	  w(28)=0.99499330568698229892987683342653327D+00
	  w(29)=0.10080994461836765285767114619375207D+01
	  w(30)=0.10212677513142154328562583032180555D+01
	  w(31)=0.10346580508699245726234039466362447D+01
	  w(32)=0.10484679850247504706572954091825522D+01
	  w(33)=0.10629558985700717776978763140505180D+01
	  w(34)=0.10784786131313242041329658604809083D+01
	  w(35)=0.10955761354344732971810572053072974D+01
	  w(36)=0.11151239584354482037298339491826482D+01
	  w(37)=0.11401180619063624455122862855205312D+01
	  w(38)=0.11835141220341260925863480224506930D+01
	endif
	if(n.eq. 39) then
	  x(1)=0.25073887894726854669302795741714363D-01
	  x(2)=0.13115344977535206560936842379305745D+00
	  x(3)=0.31833541602906467726796790884691291D+00
	  x(4)=0.58124668846864069315927281422773376D+00
	  x(5)=0.91352613422709438850688457023352385D+00
	  x(6)=0.13086737639678514000252107507549226D+01
	  x(7)=0.17605547962436960407472952283569612D+01
	  x(8)=0.22636252428551268955914110847515985D+01
	  x(9)=0.28129803717777250859910509461769834D+01
	  x(10)=0.34043140824946438449671859416412190D+01
	  x(11)=0.40338457650440728841090276546310633D+01
	  x(12)=0.46982443636267854358834483718965203D+01
	  x(13)=0.53945620989331608541306195547804236D+01
	  x(14)=0.61201809156473414574861635628622025D+01
	  x(15)=0.68727704453347486790448783722240478D+01
	  x(16)=0.76502548970767074720811251609120518D+01
	  x(17)=0.84507863323098462160487542860209942D+01
	  x(18)=0.92727223628167401159316796110942960D+01
	  x(19)=0.10114606967849358198918707785196602D+02
	  x(20)=0.10975153662712441970938925805967301D+02
	  x(21)=0.11853230621739365702183022222016007D+02
	  x(22)=0.12747847587070733510472564375959337D+02
	  x(23)=0.13658144525609888120243340381421149D+02
	  x(24)=0.14583382072903562587384840298909694D+02
	  x(25)=0.15522933859046167626161150110419840D+02
	  x(26)=0.16476280873633431411917626974172890D+02
	  x(27)=0.17443008118557223440348025178536773D+02
	  x(28)=0.18422803939836828845955096767283976D+02
	  x(29)=0.19415462654160794642166365520097315D+02
	  x(30)=0.20420891407608838363785253022797406D+02
	  x(31)=0.21439122657746818134683053358457983D+02
	  x(32)=0.22470334096155962555485530174337327D+02
	  x(33)=0.23514878337990626988585063372738659D+02
	  x(34)=0.24573326062807730352233193116262555D+02
	  x(35)=0.25646530816950477316140677430666983D+02
	  x(36)=0.26735726352941803440899093402549624D+02
	  x(37)=0.27842739408498495379262749338522553D+02
	  x(38)=0.28971263364292116904152862844057381D+02
	  x(39)=0.30124885504124957691374220303259790D+02
	  w(1)=0.64235714592269621547870883659925312D-01
	  w(2)=0.14738884967476767462457587498647626D+00
	  w(3)=0.22605882906388347075576916722638998D+00
	  w(4)=0.29868941957400840792047347349580377D+00
	  w(5)=0.36477790743588461852553450626146514D+00
	  w(6)=0.42449139500404392544297138556430582D+00
	  w(7)=0.47834607841692106111608495666587260D+00
	  w(8)=0.52697816813428965954102523028268479D+00
	  w(9)=0.57101569189262657566530378971947357D+00
	  w(10)=0.61102291998874314060685719596222043D+00
	  w(11)=0.64748637909929174227841031097341329D+00
	  w(12)=0.68082041008123117542538693669484928D+00
	  w(13)=0.71137941191494469261868971443618648D+00
	  w(14)=0.73947037235748680483737871327321045D+00
	  w(15)=0.76536316260410841572081608319422230D+00
	  w(16)=0.78929811977364650044108884685556404D+00
	  w(17)=0.81149133728208522331470931021613069D+00
	  w(18)=0.83213833324749852948798434226773679D+00
	  w(19)=0.85141671030042120893455148689099588D+00
	  w(20)=0.86948825995940692834551555279176682D+00
	  w(21)=0.88650080433177846916947828503907658D+00
	  w(22)=0.90258994713049844893504314313759096D+00
	  w(23)=0.91788083120008634452346996113192290D+00
	  w(24)=0.93248996390617033824810278019867837D+00
	  w(25)=0.94652716586289864508785285579506308D+00
	  w(26)=0.96009771620203732123144391152891330D+00
	  w(27)=0.97330480712751232807988799322629347D+00
	  w(28)=0.98625248964186451594571281020762399D+00
	  w(29)=0.99904938423317091178432747256010771D+00
	  w(30)=0.10118135423928773519719470641575754D+01
	  w(31)=0.10246790248726833638670541404280812D+01
	  w(32)=0.10378045585471802336741120598162524D+01
	  w(33)=0.10513856205994547288185003708349541D+01
	  w(34)=0.10656719105692995519518717628670856D+01
	  w(35)=0.10810033236575720128058719637920149D+01
	  w(36)=0.10978821659188529036299541985499673D+01
	  w(37)=0.11172342139069533217821117432322353D+01
	  w(38)=0.11412213696802919749018201400758699D+01
	  w(39)=0.11817543630745042104734920940245502D+01
	endif
	if(n.eq. 40) then
	  x(1)=0.24504667826819235415536368805078382D-01
	  x(2)=0.12821748571760233081917590425291564D+00
	  x(3)=0.31137268245398080912522686958254781D+00
	  x(4)=0.56890585223203204989772530097980052D+00
	  x(5)=0.89477164873848979986803442443488166D+00
	  x(6)=0.12827408059868423695348838009522296D+01
	  x(7)=0.17268891682406728449450383777730167D+01
	  x(8)=0.22218298182375098015484127245144919D+01
	  x(9)=0.27627751283186414710257849947083741D+01
	  x(10)=0.33455087714491904549163336923811585D+01
	  x(11)=0.39663215429908382780865849781548604D+01
	  x(12)=0.46219407417781592783967425930313766D+01
	  x(13)=0.53094665888658552077572494454216212D+01
	  x(14)=0.60263198727720990177658677566796541D+01
	  x(15)=0.67702004675008753409315431781578809D+01
	  x(16)=0.75390546476126569785947140189819038D+01
	  x(17)=0.83310488819235093416182280634529889D+01
	  x(18)=0.91445481836585695845087684574536979D+01
	  x(19)=0.99780976580742954951119827455841005D+01
	  x(20)=0.10830406394724086283076758263632655D+02
	  x(21)=0.11700333225833329464649068540893495D+02
	  x(22)=0.12586874116096410602949617896229029D+02
	  x(23)=0.13489151093708425577233356307260692D+02
	  x(24)=0.14406402713490455980149818060453981D+02
	  x(25)=0.15337976090804016138235965627245605D+02
	  x(26)=0.16283320582711546364862442715093493D+02
	  x(27)=0.17241983238187017946074774954468012D+02
	  x(28)=0.18213606211667723755454062484204769D+02
	  x(29)=0.19197926466344409135444948333315551D+02
	  x(30)=0.20194778344840742079213669057935476D+02
	  x(31)=0.21204100118586271861431669094599783D+02
	  x(32)=0.22225946654943975033802416874095798D+02
	  x(33)=0.23260511845374328032676203292794526D+02
	  x(34)=0.24308166596091595579309796448796988D+02
	  x(35)=0.25369517780656448024956262088380754D+02
	  x(36)=0.26445492197921929289350373437628150D+02
	  x(37)=0.27537461088154344679423957131803036D+02
	  x(38)=0.28647456688034615979177033295854926D+02
	  x(39)=0.29778334017986210824346926528960466D+02
	  x(40)=0.30938371037795256768276885850355029D+02
	  w(1)=0.62782402890552257979095429618610069D-01
	  w(2)=0.14414527700462842485507053424953483D+00
	  w(3)=0.22130189479626832804726177528209519D+00
	  w(4)=0.29274363832989835909614839692949317D+00
	  w(5)=0.35794171729349866239999755634926260D+00
	  w(6)=0.41700559913000234679714139929274097D+00
	  w(7)=0.47039246816630259084845988581946585D+00
	  w(8)=0.51869148658738262192002821393543854D+00
	  w(9)=0.56249754937039597724179884608020075D+00
	  w(10)=0.60235274236903557643074691441142932D+00
	  w(11)=0.63872841833503812836170254740864038D+00
	  w(12)=0.67202755341193376459330011130077764D+00
	  w(13)=0.70259495805650751254489705388550647D+00
	  w(14)=0.73072879442495675128554921684553847D+00
	  w(15)=0.75669053960591947660674350117915310D+00
	  w(16)=0.78071258971543677773752278881147504D+00
	  w(17)=0.80300368667774868480790928515489213D+00
	  w(18)=0.82375270171335557289893358756671660D+00
	  w(19)=0.84313133222241798403473467260482721D+00
	  w(20)=0.86129615643818069337100951088359579D+00
	  w(21)=0.87839035105925988844433049962390214D+00
	  w(22)=0.89454525993312516973787751339841634D+00
	  w(23)=0.90988192085007990606015937373740599D+00
	  w(24)=0.92451261055897526741631509139551781D+00
	  w(25)=0.93854244741608816582356666913256049D+00
	  w(26)=0.95207109077029727917107493340154178D+00
	  w(27)=0.96519459157104703361795827731839381D+00
	  w(28)=0.97800749098771688760933784578810446D+00
	  w(29)=0.99060534542288103398988141634617932D+00
	  w(30)=0.10030880346103638878219044272555038D+01
	  w(31)=0.10155646343432478406754171373904683D+01
	  w(32)=0.10281610016829541009997228684369475D+01
	  w(33)=0.10410328462714633079855275354930200D+01
	  w(34)=0.10543834318522862236733317331527360D+01
	  w(35)=0.10684925996095646638650578097440302D+01
	  w(36)=0.10837428943236429024210565330577083D+01
	  w(37)=0.11007132595541071484745998532162048D+01
	  w(38)=0.11204089142133897283315491222310811D+01
	  w(39)=0.11446152485325957748329983587609604D+01
	  w(40)=0.11821089382133427303500639027333818D+01
	endif
	if(n.eq. 41) then
	  x(1)=0.23984523421017991950954240110149840D-01
	  x(2)=0.12552961049068050991550649087002967D+00
	  x(3)=0.30497904043054141487445463098993059D+00
	  x(4)=0.55753159882308778172443908260902390D+00
	  x(5)=0.87741720264378519900105857232119888D+00
	  x(6)=0.12586491038184279833700429662712850D+01
	  x(7)=0.16954977461019282714005385059863329D+01
	  x(8)=0.21827260342350083277551675564609468D+01
	  x(9)=0.27156613618673488552701655862620100D+01
	  x(10)=0.32901776514269975848492322256788611D+01
	  x(11)=0.39026386408220425927595442772144452D+01
	  x(12)=0.45498318829865622703323424502741545D+01
	  x(13)=0.52289076040329840466824862232897431D+01
	  x(14)=0.59373274690486210758422203070949763D+01
	  x(15)=0.66728236147476049922033780603669584D+01
	  x(16)=0.74333663342336304680202374584041536D+01
	  x(17)=0.82171383220225351351473364047706127D+01
	  x(18)=0.90225136171926507699936337303370237D+01
	  x(19)=0.98480398528942885150172514840960503D+01
	  x(20)=0.10692422891002696161422136356122792D+02
	  x(21)=0.11554513294398519818173554085660726D+02
	  x(22)=0.12433294349000162171137162658851594D+02
	  x(23)=0.13327871512316947999465810426045209D+02
	  x(24)=0.14237463264860524603250269137788564D+02
	  x(25)=0.15161393403778083310839974728878587D+02
	  x(26)=0.16099084867358556039107497781515121D+02
	  x(27)=0.17050055232336600141707094735465944D+02
	  x(28)=0.18013914092609461192751041380688548D+02
	  x(29)=0.18990362612595298230644402792677283D+02
	  x(30)=0.19979195621968958818115424946881831D+02
	  x(31)=0.20980306596836737043076936970464885D+02
	  x(32)=0.21993695652042131172265726490877569D+02
	  x(33)=0.23019480134884943822726199869066477D+02
	  x(34)=0.24057907429698456525102301384322345D+02
	  x(35)=0.25109372769166562022746802540495992D+02
	  x(36)=0.26174451464588528892818430904299021D+02
	  x(37)=0.27253979420866915717169831623323262D+02
	  x(38)=0.28349199152628795417285800795070827D+02
	  x(39)=0.29461964303315493651780343498103321D+02
	  x(40)=0.30595867666822844199714381829835474D+02
	  x(41)=0.31760614049108848888636202900670469D+02
	  w(1)=0.61453777689078850632320438762690173D-01
	  w(2)=0.14116901928623146189600845445966115D+00
	  w(3)=0.21691178102955907536220081510691671D+00
	  w(4)=0.28721914904195389173224839396425523D+00
	  w(5)=0.35154646142852036261317039134155493D+00
	  w(6)=0.40995869480166080567329345285543241D+00
	  w(7)=0.46286472101062492434664363827323541D+00
	  w(8)=0.51081350768430644482265279293642379D+00
	  w(9)=0.55437048040962799699826746291364543D+00
	  w(10)=0.59405718617497671463212327580549754D+00
	  w(11)=0.63033041409925372189348991014412604D+00
	  w(12)=0.66358196163380622767391514571500011D+00
	  w(13)=0.69414710302202875524812952789943665D+00
	  w(14)=0.72231518314433129646090492315124720D+00
	  w(15)=0.74833924850753152924198730033822358D+00
	  w(16)=0.77244365386107194382248053443618119D+00
	  w(17)=0.79482962192212569529914389931946062D+00
	  w(18)=0.81567916211028024253693047285196371D+00
	  w(19)=0.83515784442209617388641618163092062D+00
	  w(20)=0.85341685831049940436088263595593162D+00
	  w(21)=0.87059467195834050468050691051757894D+00
	  w(22)=0.88681849942647339446466503432020545D+00
	  w(23)=0.90220570265884358907015894146752544D+00
	  w(24)=0.91686520534588789033136890793684870D+00
	  w(25)=0.93089897139517174817768818684271537D+00
	  w(26)=0.94440359709427723622354733379324898D+00
	  w(27)=0.95747207569971692908694649304379709D+00
	  w(28)=0.97019581022582634943773882696405053D+00
	  w(29)=0.98266696542779785605858933195122518D+00
	  w(30)=0.99498119901604997750155234825797379D+00
	  w(31)=0.10072407235150886428698413510574028D+01
	  w(32)=0.10195570153317672357218270917655900D+01
	  w(33)=0.10320537595879675230037264554994181D+01
	  w(34)=0.10448685551989436159914248491986655D+01
	  w(35)=0.10581641602642044563964418557588942D+01
	  w(36)=0.10721519908466690118586939206579700D+01
	  w(37)=0.10871725696261811400233909807866439D+01
	  w(38)=0.11037689531946959231589744376833551D+01
	  w(39)=0.11228513333944662999641650458215736D+01
	  w(40)=0.11475919063733051039832844253396615D+01
	  w(41)=0.11746124972473792613669729689718224D+01
	endif
	if(n.eq. 42) then
	  x(1)=0.23399519154861270209266876918263733D-01
	  x(2)=0.12250727755965652721048542161952355D+00
	  x(3)=0.29779385079556730531180619436781853D+00
	  x(4)=0.54476244056479006250270913369604386D+00
	  x(5)=0.85796593033227286362318864121334627D+00
	  x(6)=0.12317064415078076144993701745988801D+01
	  x(7)=0.16604889563747902148094226504326798D+01
	  x(8)=0.21392575607345967192429725400870666D+01
	  x(9)=0.26634773326195908538238654728047550D+01
	  x(10)=0.32291274032903434232366635114885867D+01
	  x(11)=0.38326531517021260597744003462139517D+01
	  x(12)=0.44709063426711255218037877057213336D+01
	  x(13)=0.51410878286817593973978546273428947D+01
	  x(14)=0.58406986322303806602462827868293971D+01
	  x(15)=0.65675004483780758945954403316136450D+01
	  x(16)=0.73194844194713235552285368612501770D+01
	  x(17)=0.80948463450017218434595633880235255D+01
	  x(18)=0.88919665623227750472779007395729423D+01
	  x(19)=0.97093931113621838591143387020565569D+01
	  x(20)=0.10545827225370498325673906947486103D+02
	  x(21)=0.11400110553885486552871952881105244D+02
	  x(22)=0.12271213789242862546302603732328862D+02
	  x(23)=0.13158226541458786229554789315443486D+02
	  x(24)=0.14060348411404792301482302718795836D+02
	  x(25)=0.14976881274013200240347032377030700D+02
	  x(26)=0.15907222823958095503371623635757715D+02
	  x(27)=0.16850861473230875020590247004292905D+02
	  x(28)=0.17807372741036434860006920644082129D+02
	  x(29)=0.18776417356958841509140256675891578D+02
	  x(30)=0.19757741446668578788603554130531847D+02
	  x(31)=0.20751179385310333458392051397822797D+02
	  x(32)=0.21756660233976980123316025128588080D+02
	  x(33)=0.22774219302597114023001267923973501D+02
	  x(34)=0.23804016061690816741247544996440411D+02
	  x(35)=0.24846361937852712742369476472958922D+02
	  x(36)=0.25901759563898949068061483558267355D+02
	  x(37)=0.26970970182886965460511419223621488D+02
	  x(38)=0.28055109231405054970309720374643803D+02
	  x(39)=0.29155906125746835044765248312614858D+02
	  x(40)=0.30276036870974266435041499789804220D+02
	  x(41)=0.31422744585242831760751869296655059D+02
	  x(42)=0.32611636111093410761441191425547004D+02
	  w(1)=0.59959549528834581166858441747535835D-01
	  w(2)=0.13782354106246585967454620913485996D+00
	  w(3)=0.21198373628310368532723373391490895D+00
	  w(4)=0.28103396142753983566464626164815854D+00
	  w(5)=0.34441587073831536613965909054968506D+00
	  w(6)=0.40214554766643612060406098862586077D+00
	  w(7)=0.45457562783648930437863100451068021D+00
	  w(8)=0.50220599007960209370082793611800298D+00
	  w(9)=0.54556427189496681240399311718647368D+00
	  w(10)=0.58514446288410426610937520308652893D+00
	  w(11)=0.62138329530779745635271638093399815D+00
	  w(12)=0.65465737881544439069614327308954671D+00
	  w(13)=0.68528975194393804137860115588409826D+00
	  w(14)=0.71355932009544631444697415645350702D+00
	  w(15)=0.73970991145305309277802052747574635D+00
	  w(16)=0.76395766265738540568008829723112285D+00
	  w(17)=0.78649652296777516902182014746358618D+00
	  w(18)=0.80750215877814202070084093065815978D+00
	  w(19)=0.82713469038170439606716399794095196D+00
	  w(20)=0.84554066844173625394631699236924760D+00
	  w(21)=0.86285460596226726259772021876415238D+00
	  w(22)=0.87920028191136079254164314988884144D+00
	  w(23)=0.89469195141805224480435754230711609D+00
	  w(24)=0.90943554210252131042580003850162029D+00
	  w(25)=0.92352988433189653338928337689139880D+00
	  w(26)=0.93706801176950760279282803821843117D+00
	  w(27)=0.95013857293707892903711353937978856D+00
	  w(28)=0.96282741662106641200580270378850400D+00
	  w(29)=0.97521945976814672363275349198374897D+00
	  w(30)=0.98740101730765483800666970637394115D+00
	  w(31)=0.99946288622088297337597850855672732D+00
	  w(32)=0.10115044760703817416924721328541636D+01
	  w(33)=0.10236398737916996282137915841303766D+01
	  w(34)=0.10360055297867141188561390663380735D+01
	  w(35)=0.10487706272210348679152502882061526D+01
	  w(36)=0.10621564297165377865894697606563568D+01
	  w(37)=0.10764703972462128511722312396159396D+01
	  w(38)=0.10921305609416507298448095752974041D+01
	  w(39)=0.11099745940549852196710389762301929D+01
	  w(40)=0.11308477326180361721696954191429541D+01
	  w(41)=0.11659271347310256050633370250579901D+01
	  w(42)=0.13123354895155796562278283090563491D+01
	endif
	end
c-----------------------------------------------    
	subroutine numthetafour(numtets,nlams)
	integer *8 nlams,numtets(nlams)
c
c This routine returns the number of Fourier modes needed in the
c phi integral for each of the discrete lambda values given
c by Norman's quadratures. (see weights.f).
c
c Input arguments:
c	nlams - number of nodes in the lambda quadrature.  This must
c	  be an integer in the range [2,39].
c
c Output arguments:
c	numtest(i) - number of Fourier modes needed for phi
c                    integral with lambda_i.
c
c Approximate accuracies of the quadrature:
c
c	NLAMS	maximum error
c
c        2     0.15318E+00
c        3     0.76505E-01
c        4     0.32149E-01
c        5     0.15630E-01
c        6     0.75110E-02
c        7     0.35030E-02
c        8     0.16243E-02
c        9     0.72230E-03
c       10     0.33074E-03
c       11     0.15035E-03
c       12     0.70952E-04
c       13     0.31751E-04
c       14     0.14589E-04
c       15     0.64300E-05
c       16     0.29477E-05
c       17     0.13222E-05
c       18     0.61488E-06
c       19     0.27435E-06
c       20     0.12534E-06
c       21     0.55324E-07
c       22     0.25257E-07
c       23     0.11293E-07
c       24     0.52063E-08
c       25     0.23256E-08
c       26     0.10580E-08
c       27     0.46835E-09
c       28     0.21286E-09
c       29     0.95164E-10
c       30     0.43599E-10
c       31     0.19516E-10
c       32     0.88491E-11
c       33     0.39313E-11
c       34     0.17821E-11
c       35     0.79603E-12
c       36     0.36460E-12
c       37     0.16331E-12
c       38     0.73497E-13
c       39     0.31530E-13
c
      if(nlams.eq. 2) then
         numtets( 1) =  3
         numtets( 2) =  1
      endif
      if(nlams.eq. 3) then
         numtets( 1) =  3
         numtets( 2) =  4
         numtets( 3) =  1
      endif
      if(nlams.eq. 4) then
         numtets( 1) =  3
         numtets( 2) =  7
         numtets( 3) =  5
         numtets( 4) =  1
      endif
      if(nlams.eq. 5) then
         numtets( 1) =  3
         numtets( 2) =  7
         numtets( 3) =  9
         numtets( 4) =  8
         numtets( 5) =  1
      endif
      if(nlams.eq. 6) then
         numtets( 1) =  3
         numtets( 2) =  7
         numtets( 3) = 10
         numtets( 4) = 12
         numtets( 5) =  8
         numtets( 6) =  1
      endif
      if(nlams.eq. 7) then
         numtets( 1) =  3
         numtets( 2) =  7
         numtets( 3) = 11
         numtets( 4) = 15
         numtets( 5) = 16
         numtets( 6) =  9
         numtets( 7) =  1
      endif
      if(nlams.eq. 8) then
         numtets( 1) =  3
         numtets( 2) =  7
         numtets( 3) = 11
         numtets( 4) = 15
         numtets( 5) = 19
         numtets( 6) = 19
         numtets( 7) =  7
         numtets( 8) =  1
      endif
      if(nlams.eq. 9) then
         numtets( 1) =  4
         numtets( 2) =  7
         numtets( 3) = 11
         numtets( 4) = 15
         numtets( 5) = 20
         numtets( 6) = 20
         numtets( 7) = 24
         numtets( 8) =  7
         numtets( 9) =  1
      endif
      if(nlams.eq.10) then
         numtets( 1) =  4
         numtets( 2) =  7
         numtets( 3) = 11
         numtets( 4) = 15
         numtets( 5) = 20
         numtets( 6) = 24
         numtets( 7) = 27
         numtets( 8) = 25
         numtets( 9) =  6
         numtets(10) =  1
      endif
      if(nlams.eq.11) then
         numtets( 1) =  4
         numtets( 2) =  7
         numtets( 3) = 11
         numtets( 4) = 16
         numtets( 5) = 20
         numtets( 6) = 24
         numtets( 7) = 28
         numtets( 8) = 32
         numtets( 9) = 29
         numtets(10) =  6
         numtets(11) =  1
      endif
      if(nlams.eq.12) then
         numtets( 1) =  4
         numtets( 2) =  8
         numtets( 3) = 11
         numtets( 4) = 16
         numtets( 5) = 20
         numtets( 6) = 24
         numtets( 7) = 28
         numtets( 8) = 33
         numtets( 9) = 35
         numtets(10) = 32
         numtets(11) =  5
         numtets(12) =  1
      endif
      if(nlams.eq.13) then
         numtets( 1) =  4
         numtets( 2) =  8
         numtets( 3) = 12
         numtets( 4) = 16
         numtets( 5) = 19
         numtets( 6) = 24
         numtets( 7) = 28
         numtets( 8) = 32
         numtets( 9) = 37
         numtets(10) = 38
         numtets(11) = 37
         numtets(12) =  6
         numtets(13) =  1
      endif
      if(nlams.eq.14) then
         numtets( 1) =  4
         numtets( 2) =  8
         numtets( 3) = 12
         numtets( 4) = 16
         numtets( 5) = 20
         numtets( 6) = 24
         numtets( 7) = 29
         numtets( 8) = 33
         numtets( 9) = 38
         numtets(10) = 40
         numtets(11) = 44
         numtets(12) = 35
         numtets(13) =  5
         numtets(14) =  1
      endif
      if(nlams.eq.15) then
         numtets( 1) =  4
         numtets( 2) =  8
         numtets( 3) = 12
         numtets( 4) = 16
         numtets( 5) = 20
         numtets( 6) = 25
         numtets( 7) = 29
         numtets( 8) = 33
         numtets( 9) = 38
         numtets(10) = 41
         numtets(11) = 45
         numtets(12) = 47
         numtets(13) = 40
         numtets(14) =  5
         numtets(15) =  1
      endif
      if(nlams.eq.16) then
         numtets( 1) =  5
         numtets( 2) =  8
         numtets( 3) = 12
         numtets( 4) = 16
         numtets( 5) = 20
         numtets( 6) = 25
         numtets( 7) = 29
         numtets( 8) = 33
         numtets( 9) = 37
         numtets(10) = 42
         numtets(11) = 47
         numtets(12) = 49
         numtets(13) = 49
         numtets(14) = 46
         numtets(15) =  6
         numtets(16) =  1
      endif
      if(nlams.eq.17) then
         numtets( 1) =  5
         numtets( 2) =  8
         numtets( 3) = 12
         numtets( 4) = 16
         numtets( 5) = 20
         numtets( 6) = 25
         numtets( 7) = 29
         numtets( 8) = 34
         numtets( 9) = 38
         numtets(10) = 43
         numtets(11) = 47
         numtets(12) = 51
         numtets(13) = 56
         numtets(14) = 56
         numtets(15) = 49
         numtets(16) =  5
         numtets(17) =  1
      endif
      if(nlams.eq.18) then
         numtets( 1) =  5
         numtets( 2) =  8
         numtets( 3) = 12
         numtets( 4) = 16
         numtets( 5) = 20
         numtets( 6) = 25
         numtets( 7) = 29
         numtets( 8) = 34
         numtets( 9) = 38
         numtets(10) = 43
         numtets(11) = 47
         numtets(12) = 51
         numtets(13) = 56
         numtets(14) = 59
         numtets(15) = 59
         numtets(16) = 51
         numtets(17) =  4
         numtets(18) =  1
      endif
      if(nlams.eq.19) then
         numtets( 1) =  5
         numtets( 2) =  9
         numtets( 3) = 12
         numtets( 4) = 16
         numtets( 5) = 21
         numtets( 6) = 25
         numtets( 7) = 29
         numtets( 8) = 34
         numtets( 9) = 38
         numtets(10) = 43
         numtets(11) = 47
         numtets(12) = 52
         numtets(13) = 56
         numtets(14) = 60
         numtets(15) = 65
         numtets(16) = 65
         numtets(17) = 53
         numtets(18) =  5
         numtets(19) =  1
      endif
      if(nlams.eq.20) then
         numtets( 1) =  5
         numtets( 2) =  9
         numtets( 3) = 13
         numtets( 4) = 17
         numtets( 5) = 21
         numtets( 6) = 25
         numtets( 7) = 29
         numtets( 8) = 34
         numtets( 9) = 38
         numtets(10) = 43
         numtets(11) = 47
         numtets(12) = 52
         numtets(13) = 56
         numtets(14) = 60
         numtets(15) = 64
         numtets(16) = 68
         numtets(17) = 66
         numtets(18) = 53
         numtets(19) =  4
         numtets(20) =  1
      endif
      if(nlams.eq.21) then
         numtets( 1) =  5
         numtets( 2) =  9
         numtets( 3) = 13
         numtets( 4) = 17
         numtets( 5) = 21
         numtets( 6) = 25
         numtets( 7) = 30
         numtets( 8) = 34
         numtets( 9) = 38
         numtets(10) = 43
         numtets(11) = 48
         numtets(12) = 52
         numtets(13) = 57
         numtets(14) = 61
         numtets(15) = 65
         numtets(16) = 68
         numtets(17) = 71
         numtets(18) = 68
         numtets(19) = 54
         numtets(20) =  5
         numtets(21) =  1
      endif
      if(nlams.eq.22) then
         numtets( 1) =  6
         numtets( 2) =  9
         numtets( 3) = 13
         numtets( 4) = 17
         numtets( 5) = 21
         numtets( 6) = 25
         numtets( 7) = 30
         numtets( 8) = 34
         numtets( 9) = 39
         numtets(10) = 43
         numtets(11) = 48
         numtets(12) = 52
         numtets(13) = 57
         numtets(14) = 61
         numtets(15) = 66
         numtets(16) = 68
         numtets(17) = 74
         numtets(18) = 74
         numtets(19) = 75
         numtets(20) = 63
         numtets(21) =  4
         numtets(22) =  1
      endif
      if(nlams.eq.23) then
         numtets( 1) =  6
         numtets( 2) =  9
         numtets( 3) = 13
         numtets( 4) = 17
         numtets( 5) = 21
         numtets( 6) = 25
         numtets( 7) = 30
         numtets( 8) = 34
         numtets( 9) = 39
         numtets(10) = 43
         numtets(11) = 48
         numtets(12) = 52
         numtets(13) = 57
         numtets(14) = 61
         numtets(15) = 66
         numtets(16) = 70
         numtets(17) = 75
         numtets(18) = 78
         numtets(19) = 82
         numtets(20) = 78
         numtets(21) = 67
         numtets(22) =  4
         numtets(23) =  1
      endif
      if(nlams.eq.24) then
         numtets( 1) =  6
         numtets( 2) =  9
         numtets( 3) = 13
         numtets( 4) = 17
         numtets( 5) = 21
         numtets( 6) = 25
         numtets( 7) = 30
         numtets( 8) = 34
         numtets( 9) = 39
         numtets(10) = 43
         numtets(11) = 48
         numtets(12) = 52
         numtets(13) = 57
         numtets(14) = 61
         numtets(15) = 66
         numtets(16) = 70
         numtets(17) = 75
         numtets(18) = 79
         numtets(19) = 82
         numtets(20) = 83
         numtets(21) = 83
         numtets(22) = 67
         numtets(23) =  5
         numtets(24) =  1
      endif
      if(nlams.eq.25) then
         numtets( 1) =  6
         numtets( 2) = 10
         numtets( 3) = 13
         numtets( 4) = 17
         numtets( 5) = 21
         numtets( 6) = 25
         numtets( 7) = 30
         numtets( 8) = 34
         numtets( 9) = 39
         numtets(10) = 43
         numtets(11) = 48
         numtets(12) = 52
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 66
         numtets(16) = 71
         numtets(17) = 75
         numtets(18) = 79
         numtets(19) = 83
         numtets(20) = 88
         numtets(21) = 88
         numtets(22) = 89
         numtets(23) = 71
         numtets(24) =  4
         numtets(25) =  1
      endif
      if(nlams.eq.26) then
         numtets( 1) =  6
         numtets( 2) = 10
         numtets( 3) = 13
         numtets( 4) = 17
         numtets( 5) = 21
         numtets( 6) = 26
         numtets( 7) = 30
         numtets( 8) = 34
         numtets( 9) = 39
         numtets(10) = 43
         numtets(11) = 48
         numtets(12) = 52
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 66
         numtets(16) = 71
         numtets(17) = 75
         numtets(18) = 80
         numtets(19) = 84
         numtets(20) = 88
         numtets(21) = 91
         numtets(22) = 92
         numtets(23) = 89
         numtets(24) = 69
         numtets(25) =  5
         numtets(26) =  1
      endif
      if(nlams.eq.27) then
         numtets( 1) =  6
         numtets( 2) = 10
         numtets( 3) = 14
         numtets( 4) = 17
         numtets( 5) = 22
         numtets( 6) = 26
         numtets( 7) = 30
         numtets( 8) = 34
         numtets( 9) = 39
         numtets(10) = 43
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 66
         numtets(16) = 71
         numtets(17) = 75
         numtets(18) = 80
         numtets(19) = 84
         numtets(20) = 89
         numtets(21) = 92
         numtets(22) = 95
         numtets(23) = 97
         numtets(24) = 95
         numtets(25) = 75
         numtets(26) =  4
         numtets(27) =  1
      endif
      if(nlams.eq.28) then
         numtets( 1) =  6
         numtets( 2) = 10
         numtets( 3) = 14
         numtets( 4) = 18
         numtets( 5) = 22
         numtets( 6) = 26
         numtets( 7) = 30
         numtets( 8) = 34
         numtets( 9) = 39
         numtets(10) = 43
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 66
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 80
         numtets(19) = 85
         numtets(20) = 89
         numtets(21) = 93
         numtets(22) = 98
         numtets(23) =102
         numtets(24) =102
         numtets(25) = 99
         numtets(26) = 80
         numtets(27) =  4
         numtets(28) =  1
      endif
      if(nlams.eq.29) then
         numtets( 1) =  6
         numtets( 2) = 10
         numtets( 3) = 14
         numtets( 4) = 18
         numtets( 5) = 22
         numtets( 6) = 26
         numtets( 7) = 30
         numtets( 8) = 35
         numtets( 9) = 39
         numtets(10) = 43
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 66
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 80
         numtets(19) = 85
         numtets(20) = 89
         numtets(21) = 94
         numtets(22) = 98
         numtets(23) =102
         numtets(24) =103
         numtets(25) =108
         numtets(26) =101
         numtets(27) = 86
         numtets(28) =  4
         numtets(29) =  1
      endif
      if(nlams.eq.30) then
         numtets( 1) =  7
         numtets( 2) = 10
         numtets( 3) = 14
         numtets( 4) = 18
         numtets( 5) = 22
         numtets( 6) = 26
         numtets( 7) = 30
         numtets( 8) = 35
         numtets( 9) = 39
         numtets(10) = 43
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 66
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 80
         numtets(19) = 85
         numtets(20) = 89
         numtets(21) = 94
         numtets(22) = 97
         numtets(23) =103
         numtets(24) =107
         numtets(25) =110
         numtets(26) =112
         numtets(27) =108
         numtets(28) = 84
         numtets(29) =  4
         numtets(30) =  1
      endif
      if(nlams.eq.31) then
         numtets( 1) =  7
         numtets( 2) = 10
         numtets( 3) = 14
         numtets( 4) = 18
         numtets( 5) = 22
         numtets( 6) = 26
         numtets( 7) = 30
         numtets( 8) = 35
         numtets( 9) = 39
         numtets(10) = 44
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 67
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 80
         numtets(19) = 85
         numtets(20) = 90
         numtets(21) = 94
         numtets(22) = 99
         numtets(23) =103
         numtets(24) =107
         numtets(25) =110
         numtets(26) =115
         numtets(27) =117
         numtets(28) =111
         numtets(29) = 85
         numtets(30) =  5
         numtets(31) =  1
      endif
      if(nlams.eq.32) then
         numtets( 1) =  7
         numtets( 2) = 11
         numtets( 3) = 14
         numtets( 4) = 18
         numtets( 5) = 22
         numtets( 6) = 26
         numtets( 7) = 30
         numtets( 8) = 35
         numtets( 9) = 39
         numtets(10) = 44
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 67
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 80
         numtets(19) = 85
         numtets(20) = 90
         numtets(21) = 94
         numtets(22) = 99
         numtets(23) =103
         numtets(24) =107
         numtets(25) =112
         numtets(26) =114
         numtets(27) =119
         numtets(28) =117
         numtets(29) =108
         numtets(30) = 89
         numtets(31) =  4
         numtets(32) =  1
      endif
      if(nlams.eq.33) then
         numtets( 1) =  7
         numtets( 2) = 11
         numtets( 3) = 15
         numtets( 4) = 18
         numtets( 5) = 22
         numtets( 6) = 26
         numtets( 7) = 31
         numtets( 8) = 35
         numtets( 9) = 39
         numtets(10) = 44
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 67
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 81
         numtets(19) = 85
         numtets(20) = 90
         numtets(21) = 94
         numtets(22) = 99
         numtets(23) =103
         numtets(24) =107
         numtets(25) =112
         numtets(26) =116
         numtets(27) =120
         numtets(28) =124
         numtets(29) =124
         numtets(30) =116
         numtets(31) = 93
         numtets(32) =  5
         numtets(33) =  1
      endif
      if(nlams.eq.34) then
         numtets( 1) =  7
         numtets( 2) = 11
         numtets( 3) = 15
         numtets( 4) = 18
         numtets( 5) = 22
         numtets( 6) = 26
         numtets( 7) = 31
         numtets( 8) = 35
         numtets( 9) = 39
         numtets(10) = 44
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 67
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 81
         numtets(19) = 85
         numtets(20) = 90
         numtets(21) = 94
         numtets(22) = 99
         numtets(23) =104
         numtets(24) =108
         numtets(25) =112
         numtets(26) =117
         numtets(27) =121
         numtets(28) =125
         numtets(29) =128
         numtets(30) =131
         numtets(31) =127
         numtets(32) = 95
         numtets(33) =  4
         numtets(34) =  1
      endif
      if(nlams.eq.35) then
         numtets( 1) =  7
         numtets( 2) = 11
         numtets( 3) = 15
         numtets( 4) = 19
         numtets( 5) = 23
         numtets( 6) = 27
         numtets( 7) = 31
         numtets( 8) = 35
         numtets( 9) = 39
         numtets(10) = 44
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 67
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 81
         numtets(19) = 85
         numtets(20) = 90
         numtets(21) = 95
         numtets(22) = 99
         numtets(23) =104
         numtets(24) =108
         numtets(25) =113
         numtets(26) =117
         numtets(27) =121
         numtets(28) =125
         numtets(29) =130
         numtets(30) =133
         numtets(31) =133
         numtets(32) =124
         numtets(33) =105
         numtets(34) =  4
         numtets(35) =  1
      endif
      if(nlams.eq.36) then
         numtets( 1) =  7
         numtets( 2) = 11
         numtets( 3) = 15
         numtets( 4) = 19
         numtets( 5) = 23
         numtets( 6) = 27
         numtets( 7) = 31
         numtets( 8) = 35
         numtets( 9) = 39
         numtets(10) = 44
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 67
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 81
         numtets(19) = 85
         numtets(20) = 90
         numtets(21) = 95
         numtets(22) = 99
         numtets(23) =104
         numtets(24) =108
         numtets(25) =113
         numtets(26) =117
         numtets(27) =122
         numtets(28) =126
         numtets(29) =130
         numtets(30) =133
         numtets(31) =137
         numtets(32) =137
         numtets(33) =128
         numtets(34) =103
         numtets(35) =  4
         numtets(36) =  1
      endif
      if(nlams.eq.37) then
         numtets( 1) =  8
         numtets( 2) = 11
         numtets( 3) = 15
         numtets( 4) = 19
         numtets( 5) = 23
         numtets( 6) = 27
         numtets( 7) = 31
         numtets( 8) = 35
         numtets( 9) = 40
         numtets(10) = 44
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 67
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 81
         numtets(19) = 85
         numtets(20) = 90
         numtets(21) = 95
         numtets(22) = 99
         numtets(23) =104
         numtets(24) =108
         numtets(25) =113
         numtets(26) =118
         numtets(27) =122
         numtets(28) =126
         numtets(29) =131
         numtets(30) =135
         numtets(31) =139
         numtets(32) =142
         numtets(33) =144
         numtets(34) =139
         numtets(35) =107
         numtets(36) =  4
         numtets(37) =  1
      endif
      if(nlams.eq.38) then
         numtets( 1) =  8
         numtets( 2) = 12
         numtets( 3) = 15
         numtets( 4) = 19
         numtets( 5) = 23
         numtets( 6) = 27
         numtets( 7) = 31
         numtets( 8) = 35
         numtets( 9) = 40
         numtets(10) = 44
         numtets(11) = 48
         numtets(12) = 53
         numtets(13) = 57
         numtets(14) = 62
         numtets(15) = 67
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 81
         numtets(19) = 85
         numtets(20) = 90
         numtets(21) = 95
         numtets(22) = 99
         numtets(23) =104
         numtets(24) =109
         numtets(25) =113
         numtets(26) =118
         numtets(27) =122
         numtets(28) =127
         numtets(29) =131
         numtets(30) =135
         numtets(31) =139
         numtets(32) =143
         numtets(33) =143
         numtets(34) =141
         numtets(35) =139
         numtets(36) =106
         numtets(37) =  4
         numtets(38) =  1
      endif
      if(nlams.eq.39) then
         numtets( 1) =  8
         numtets( 2) = 12
         numtets( 3) = 15
         numtets( 4) = 19
         numtets( 5) = 23
         numtets( 6) = 27
         numtets( 7) = 31
         numtets( 8) = 35
         numtets( 9) = 40
         numtets(10) = 44
         numtets(11) = 49
         numtets(12) = 53
         numtets(13) = 58
         numtets(14) = 62
         numtets(15) = 67
         numtets(16) = 71
         numtets(17) = 76
         numtets(18) = 81
         numtets(19) = 85
         numtets(20) = 90
         numtets(21) = 95
         numtets(22) = 99
         numtets(23) =104
         numtets(24) =109
         numtets(25) =113
         numtets(26) =118
         numtets(27) =122
         numtets(28) =127
         numtets(29) =131
         numtets(30) =136
         numtets(31) =140
         numtets(32) =142
         numtets(33) =147
         numtets(34) =151
         numtets(35) =151
         numtets(36) =139
         numtets(37) =111
         numtets(38) =  5
         numtets(39) =  1
      endif
      return
      end
c----------------------------------      
	subroutine numthetahalf(numtets,nlams)
	integer *8 nlams,numtets(nlams)
c
c This routine returns the number of Fourier modes needed in the
c phi integral for each of the discrete lambda values given
c by Norman's quadratures. (see weights.f).
c
c Input arguments:
c	nlams - number of nodes in the lambda quadrature.  This must
c	  be an integer in the range [2,39].
c
c Output arguments:
c	numtest(i) - number of Fourier modes needed for phi
c                    integral with lambda_i.
c
c Approximate accuracies of the quadrature:
c
c	NLAMS	maximum error
c
c        2     0.15318E+00
c        3     0.76505E-01
c        4     0.32149E-01
c        5     0.15630E-01
c        6     0.75110E-02
c        7     0.35030E-02
c        8     0.16243E-02
c        9     0.72230E-03
c       10     0.33074E-03
c       11     0.15035E-03
c       12     0.70952E-04
c       13     0.31751E-04
c       14     0.14589E-04
c       15     0.64300E-05
c       16     0.29477E-05
c       17     0.13222E-05
c       18     0.61488E-06
c       19     0.27435E-06
c       20     0.12534E-06
c       21     0.55324E-07
c       22     0.25257E-07
c       23     0.11293E-07
c       24     0.52063E-08
c       25     0.23256E-08
c       26     0.10580E-08
c       27     0.46835E-09
c       28     0.21286E-09
c       29     0.95164E-10
c       30     0.43599E-10
c       31     0.19516E-10
c       32     0.88491E-11
c       33     0.39313E-11
c       34     0.17821E-11
c       35     0.79603E-12
c       36     0.36460E-12
c       37     0.16331E-12
c       38     0.73497E-13
c       39     0.31530E-13
c
      if(nlams.eq. 2) then
         numtets( 1) = 1
         numtets( 2) = 1
      endif
      if(nlams.eq. 3) then
         numtets( 1) = 1
         numtets( 2) = 2
         numtets( 3) = 2
      endif
      if(nlams.eq. 4) then
         numtets( 1) = 1
         numtets( 2) = 2
         numtets( 3) = 2
         numtets( 4) = 1
      endif
      if(nlams.eq. 5) then
         numtets( 1) = 2
         numtets( 2) = 2
         numtets( 3) = 3
         numtets( 4) = 3
         numtets( 5) = 1
      endif
      if(nlams.eq. 6) then
         numtets( 1) = 2
         numtets( 2) = 3
         numtets( 3) = 3
         numtets( 4) = 3
         numtets( 5) = 2
         numtets( 6) = 1
      endif
      if(nlams.eq. 7) then
         numtets( 1) = 2
         numtets( 2) = 3
         numtets( 3) = 3
         numtets( 4) = 4
         numtets( 5) = 4
         numtets( 6) = 3
         numtets( 7) = 1
      endif
      if(nlams.eq. 8) then
         numtets( 1) = 2
         numtets( 2) = 3
         numtets( 3) = 4
         numtets( 4) = 4
         numtets( 5) = 5
         numtets( 6) = 5
         numtets( 7) = 5
         numtets( 8) = 1
      endif
      if(nlams.eq. 9) then
         numtets( 1) = 2
         numtets( 2) = 3
         numtets( 3) = 4
         numtets( 4) = 5
         numtets( 5) = 5
         numtets( 6) = 4
         numtets( 7) = 5
         numtets( 8) = 4
         numtets( 9) = 1
      endif
      if(nlams.eq.10) then
         numtets( 1) = 2
         numtets( 2) = 3
         numtets( 3) = 4
         numtets( 4) = 5
         numtets( 5) = 5
         numtets( 6) = 6
         numtets( 7) = 6
         numtets( 8) = 6
         numtets( 9) = 6
         numtets(10) = 1
      endif
      if(nlams.eq.11) then
         numtets( 1) = 2
         numtets( 2) = 3
         numtets( 3) = 4
         numtets( 4) = 5
         numtets( 5) = 6
         numtets( 6) = 6
         numtets( 7) = 6
         numtets( 8) = 7
         numtets( 9) = 7
         numtets(10) = 5
         numtets(11) = 1
      endif
      if(nlams.eq.12) then
         numtets( 1) = 2
         numtets( 2) = 4
         numtets( 3) = 5
         numtets( 4) = 5
         numtets( 5) = 6
         numtets( 6) = 6
         numtets( 7) = 7
         numtets( 8) = 7
         numtets( 9) = 7
         numtets(10) = 7
         numtets(11) = 5
         numtets(12) = 1
      endif
      if(nlams.eq.13) then
         numtets( 1) = 3
         numtets( 2) = 4
         numtets( 3) = 5
         numtets( 4) = 6
         numtets( 5) = 6
         numtets( 6) = 7
         numtets( 7) = 7
         numtets( 8) = 8
         numtets( 9) = 8
         numtets(10) = 8
         numtets(11) = 7
         numtets(12) = 5
         numtets(13) = 1
      endif
      if(nlams.eq.14) then
         numtets( 1) = 3
         numtets( 2) = 4
         numtets( 3) = 5
         numtets( 4) = 6
         numtets( 5) = 7
         numtets( 6) = 7
         numtets( 7) = 8
         numtets( 8) = 8
         numtets( 9) = 8
         numtets(10) = 9
         numtets(11) = 8
         numtets(12) = 8
         numtets(13) = 7
         numtets(14) = 1
      endif
      if(nlams.eq.15) then
         numtets( 1) = 3
         numtets( 2) = 4
         numtets( 3) = 5
         numtets( 4) = 6
         numtets( 5) = 7
         numtets( 6) = 7
         numtets( 7) = 8
         numtets( 8) = 8
         numtets( 9) = 9
         numtets(10) = 9
         numtets(11) = 9
         numtets(12) = 9
         numtets(13) = 8
         numtets(14) = 7
         numtets(15) = 1
      endif
      if(nlams.eq.16) then
         numtets( 1) = 3
         numtets( 2) = 4
         numtets( 3) = 5
         numtets( 4) = 6
         numtets( 5) = 7
         numtets( 6) = 8
         numtets( 7) = 8
         numtets( 8) = 9
         numtets( 9) = 9
         numtets(10) =10
         numtets(11) =10
         numtets(12) =10
         numtets(13) =10
         numtets(14) = 9
         numtets(15) = 8
         numtets(16) = 1
      endif
      if(nlams.eq.17) then
         numtets( 1) = 3
         numtets( 2) = 4
         numtets( 3) = 6
         numtets( 4) = 6
         numtets( 5) = 7
         numtets( 6) = 8
         numtets( 7) = 9
         numtets( 8) = 9
         numtets( 9) =10
         numtets(10) =10
         numtets(11) =10
         numtets(12) =10
         numtets(13) =11
         numtets(14) =10
         numtets(15) = 9
         numtets(16) = 8
         numtets(17) = 1
      endif
      if(nlams.eq.18) then
         numtets( 1) = 3
         numtets( 2) = 5
         numtets( 3) = 6
         numtets( 4) = 7
         numtets( 5) = 7
         numtets( 6) = 8
         numtets( 7) = 9
         numtets( 8) = 9
         numtets( 9) =10
         numtets(10) =10
         numtets(11) =11
         numtets(12) =11
         numtets(13) =11
         numtets(14) =11
         numtets(15) =11
         numtets(16) =11
         numtets(17) = 7
         numtets(18) = 1
      endif
      if(nlams.eq.19) then
         numtets( 1) = 3
         numtets( 2) = 5
         numtets( 3) = 6
         numtets( 4) = 7
         numtets( 5) = 8
         numtets( 6) = 8
         numtets( 7) = 9
         numtets( 8) =10
         numtets( 9) =10
         numtets(10) =11
         numtets(11) =11
         numtets(12) =12
         numtets(13) =12
         numtets(14) =12
         numtets(15) =12
         numtets(16) =12
         numtets(17) = 9
         numtets(18) = 9
         numtets(19) = 1
      endif
      if(nlams.eq.20) then
         numtets( 1) = 3
         numtets( 2) = 5
         numtets( 3) = 6
         numtets( 4) = 7
         numtets( 5) = 8
         numtets( 6) = 9
         numtets( 7) = 9
         numtets( 8) =10
         numtets( 9) =11
         numtets(10) =11
         numtets(11) =12
         numtets(12) =12
         numtets(13) =12
         numtets(14) =12
         numtets(15) =13
         numtets(16) =13
         numtets(17) =12
         numtets(18) =12
         numtets(19) = 8
         numtets(20) = 1
      endif
      if(nlams.eq.21) then
         numtets( 1) = 4
         numtets( 2) = 5
         numtets( 3) = 6
         numtets( 4) = 7
         numtets( 5) = 8
         numtets( 6) = 9
         numtets( 7) =10
         numtets( 8) =10
         numtets( 9) =11
         numtets(10) =12
         numtets(11) =12
         numtets(12) =12
         numtets(13) =13
         numtets(14) =13
         numtets(15) =13
         numtets(16) =13
         numtets(17) =13
         numtets(18) =13
         numtets(19) =12
         numtets(20) =10
         numtets(21) = 1
      endif
      if(nlams.eq.22) then
         numtets( 1) = 4
         numtets( 2) = 5
         numtets( 3) = 6
         numtets( 4) = 7
         numtets( 5) = 8
         numtets( 6) = 9
         numtets( 7) =10
         numtets( 8) =11
         numtets( 9) =11
         numtets(10) =12
         numtets(11) =12
         numtets(12) =13
         numtets(13) =13
         numtets(14) =13
         numtets(15) =14
         numtets(16) =14
         numtets(17) =14
         numtets(18) =14
         numtets(19) =14
         numtets(20) =12
         numtets(21) =10
         numtets(22) = 1
      endif
      if(nlams.eq.23) then
         numtets( 1) = 4
         numtets( 2) = 5
         numtets( 3) = 6
         numtets( 4) = 8
         numtets( 5) = 9
         numtets( 6) = 9
         numtets( 7) =10
         numtets( 8) =11
         numtets( 9) =12
         numtets(10) =12
         numtets(11) =13
         numtets(12) =13
         numtets(13) =14
         numtets(14) =14
         numtets(15) =14
         numtets(16) =14
         numtets(17) =15
         numtets(18) =15
         numtets(19) =14
         numtets(20) =14
         numtets(21) =14
         numtets(22) = 9
         numtets(23) = 1
      endif
      if(nlams.eq.24) then
         numtets( 1) = 4
         numtets( 2) = 5
         numtets( 3) = 7
         numtets( 4) = 8
         numtets( 5) = 9
         numtets( 6) =10
         numtets( 7) =10
         numtets( 8) =11
         numtets( 9) =12
         numtets(10) =13
         numtets(11) =13
         numtets(12) =14
         numtets(13) =14
         numtets(14) =14
         numtets(15) =15
         numtets(16) =15
         numtets(17) =15
         numtets(18) =15
         numtets(19) =15
         numtets(20) =15
         numtets(21) =15
         numtets(22) =14
         numtets(23) =10
         numtets(24) = 1
      endif
      if(nlams.eq.25) then
         numtets( 1) = 4
         numtets( 2) = 6
         numtets( 3) = 7
         numtets( 4) = 8
         numtets( 5) = 9
         numtets( 6) =10
         numtets( 7) =11
         numtets( 8) =11
         numtets( 9) =12
         numtets(10) =13
         numtets(11) =13
         numtets(12) =14
         numtets(13) =14
         numtets(14) =15
         numtets(15) =15
         numtets(16) =15
         numtets(17) =16
         numtets(18) =16
         numtets(19) =16
         numtets(20) =16
         numtets(21) =16
         numtets(22) =13
         numtets(23) =14
         numtets(24) =11
         numtets(25) = 1
      endif
      if(nlams.eq.26) then
         numtets( 1) = 4
         numtets( 2) = 6
         numtets( 3) = 7
         numtets( 4) = 8
         numtets( 5) = 9
         numtets( 6) =10
         numtets( 7) =11
         numtets( 8) =12
         numtets( 9) =12
         numtets(10) =13
         numtets(11) =14
         numtets(12) =14
         numtets(13) =15
         numtets(14) =15
         numtets(15) =16
         numtets(16) =16
         numtets(17) =16
         numtets(18) =16
         numtets(19) =17
         numtets(20) =17
         numtets(21) =17
         numtets(22) =16
         numtets(23) =15
         numtets(24) =14
         numtets(25) =11
         numtets(26) = 1
      endif
      if(nlams.eq.27) then
         numtets( 1) = 4
         numtets( 2) = 6
         numtets( 3) = 7
         numtets( 4) = 8
         numtets( 5) = 9
         numtets( 6) =10
         numtets( 7) =11
         numtets( 8) =12
         numtets( 9) =13
         numtets(10) =13
         numtets(11) =14
         numtets(12) =15
         numtets(13) =15
         numtets(14) =16
         numtets(15) =16
         numtets(16) =16
         numtets(17) =17
         numtets(18) =17
         numtets(19) =17
         numtets(20) =17
         numtets(21) =17
         numtets(22) =17
         numtets(23) =17
         numtets(24) =16
         numtets(25) =15
         numtets(26) =12
         numtets(27) = 1
      endif
      if(nlams.eq.28) then
         numtets( 1) = 4
         numtets( 2) = 6
         numtets( 3) = 7
         numtets( 4) = 8
         numtets( 5) = 9
         numtets( 6) =10
         numtets( 7) =11
         numtets( 8) =12
         numtets( 9) =13
         numtets(10) =14
         numtets(11) =14
         numtets(12) =15
         numtets(13) =15
         numtets(14) =16
         numtets(15) =16
         numtets(16) =17
         numtets(17) =17
         numtets(18) =17
         numtets(19) =18
         numtets(20) =18
         numtets(21) =18
         numtets(22) =18
         numtets(23) =18
         numtets(24) =18
         numtets(25) =17
         numtets(26) =15
         numtets(27) =12
         numtets(28) = 1
      endif
      if(nlams.eq.29) then
         numtets( 1) = 4
         numtets( 2) = 6
         numtets( 3) = 7
         numtets( 4) = 9
         numtets( 5) =10
         numtets( 6) =11
         numtets( 7) =12
         numtets( 8) =12
         numtets( 9) =13
         numtets(10) =14
         numtets(11) =15
         numtets(12) =15
         numtets(13) =16
         numtets(14) =16
         numtets(15) =17
         numtets(16) =17
         numtets(17) =18
         numtets(18) =18
         numtets(19) =18
         numtets(20) =18
         numtets(21) =19
         numtets(22) =19
         numtets(23) =19
         numtets(24) =19
         numtets(25) =18
         numtets(26) =16
         numtets(27) =17
         numtets(28) =11
         numtets(29) = 1
      endif
      if(nlams.eq.30) then
         numtets( 1) = 5
         numtets( 2) = 6
         numtets( 3) = 8
         numtets( 4) = 9
         numtets( 5) =10
         numtets( 6) =11
         numtets( 7) =12
         numtets( 8) =13
         numtets( 9) =13
         numtets(10) =14
         numtets(11) =15
         numtets(12) =16
         numtets(13) =16
         numtets(14) =17
         numtets(15) =17
         numtets(16) =18
         numtets(17) =18
         numtets(18) =18
         numtets(19) =19
         numtets(20) =19
         numtets(21) =19
         numtets(22) =19
         numtets(23) =19
         numtets(24) =19
         numtets(25) =19
         numtets(26) =19
         numtets(27) =19
         numtets(28) =15
         numtets(29) =12
         numtets(30) = 1
      endif
      if(nlams.eq.31) then
         numtets( 1) = 5
         numtets( 2) = 6
         numtets( 3) = 8
         numtets( 4) = 9
         numtets( 5) =10
         numtets( 6) =11
         numtets( 7) =12
         numtets( 8) =13
         numtets( 9) =14
         numtets(10) =14
         numtets(11) =15
         numtets(12) =16
         numtets(13) =16
         numtets(14) =17
         numtets(15) =18
         numtets(16) =18
         numtets(17) =18
         numtets(18) =19
         numtets(19) =19
         numtets(20) =19
         numtets(21) =20
         numtets(22) =20
         numtets(23) =20
         numtets(24) =20
         numtets(25) =20
         numtets(26) =20
         numtets(27) =20
         numtets(28) =18
         numtets(29) =19
         numtets(30) =11
         numtets(31) = 1
      endif
      if(nlams.eq.32) then
         numtets( 1) = 5
         numtets( 2) = 6
         numtets( 3) = 8
         numtets( 4) = 9
         numtets( 5) =10
         numtets( 6) =11
         numtets( 7) =12
         numtets( 8) =13
         numtets( 9) =14
         numtets(10) =15
         numtets(11) =15
         numtets(12) =16
         numtets(13) =17
         numtets(14) =17
         numtets(15) =18
         numtets(16) =18
         numtets(17) =19
         numtets(18) =19
         numtets(19) =20
         numtets(20) =20
         numtets(21) =20
         numtets(22) =20
         numtets(23) =21
         numtets(24) =21
         numtets(25) =21
         numtets(26) =21
         numtets(27) =20
         numtets(28) =20
         numtets(29) =19
         numtets(30) =17
         numtets(31) =13
         numtets(32) = 1
      endif
      if(nlams.eq.33) then
         numtets( 1) = 5
         numtets( 2) = 7
         numtets( 3) = 8
         numtets( 4) = 9
         numtets( 5) =10
         numtets( 6) =11
         numtets( 7) =12
         numtets( 8) =13
         numtets( 9) =14
         numtets(10) =15
         numtets(11) =16
         numtets(12) =16
         numtets(13) =17
         numtets(14) =18
         numtets(15) =18
         numtets(16) =19
         numtets(17) =19
         numtets(18) =20
         numtets(19) =20
         numtets(20) =20
         numtets(21) =21
         numtets(22) =21
         numtets(23) =21
         numtets(24) =21
         numtets(25) =21
         numtets(26) =21
         numtets(27) =20
         numtets(28) =21
         numtets(29) =19
         numtets(30) =19
         numtets(31) =18
         numtets(32) =14
         numtets(33) = 1
      endif
      if(nlams.eq.34) then
         numtets( 1) = 5
         numtets( 2) = 7
         numtets( 3) = 8
         numtets( 4) = 9
         numtets( 5) =10
         numtets( 6) =12
         numtets( 7) =13
         numtets( 8) =13
         numtets( 9) =14
         numtets(10) =15
         numtets(11) =16
         numtets(12) =17
         numtets(13) =17
         numtets(14) =18
         numtets(15) =19
         numtets(16) =19
         numtets(17) =20
         numtets(18) =20
         numtets(19) =20
         numtets(20) =21
         numtets(21) =21
         numtets(22) =21
         numtets(23) =22
         numtets(24) =22
         numtets(25) =22
         numtets(26) =22
         numtets(27) =22
         numtets(28) =22
         numtets(29) =22
         numtets(30) =21
         numtets(31) =21
         numtets(32) =18
         numtets(33) =13
         numtets(34) = 1
      endif
      if(nlams.eq.35) then
         numtets( 1) = 5
         numtets( 2) = 7
         numtets( 3) = 8
         numtets( 4) = 9
         numtets( 5) =11
         numtets( 6) =12
         numtets( 7) =13
         numtets( 8) =14
         numtets( 9) =15
         numtets(10) =15
         numtets(11) =16
         numtets(12) =17
         numtets(13) =18
         numtets(14) =18
         numtets(15) =19
         numtets(16) =19
         numtets(17) =20
         numtets(18) =20
         numtets(19) =21
         numtets(20) =21
         numtets(21) =22
         numtets(22) =22
         numtets(23) =22
         numtets(24) =22
         numtets(25) =23
         numtets(26) =23
         numtets(27) =23
         numtets(28) =23
         numtets(29) =23
         numtets(30) =22
         numtets(31) =22
         numtets(32) =21
         numtets(33) =19
         numtets(34) =13
         numtets(35) = 1
      endif
      if(nlams.eq.36) then
         numtets( 1) = 5
         numtets( 2) = 7
         numtets( 3) = 8
         numtets( 4) =10
         numtets( 5) =11
         numtets( 6) =12
         numtets( 7) =13
         numtets( 8) =14
         numtets( 9) =15
         numtets(10) =16
         numtets(11) =16
         numtets(12) =17
         numtets(13) =18
         numtets(14) =19
         numtets(15) =19
         numtets(16) =20
         numtets(17) =20
         numtets(18) =21
         numtets(19) =21
         numtets(20) =22
         numtets(21) =22
         numtets(22) =22
         numtets(23) =23
         numtets(24) =23
         numtets(25) =23
         numtets(26) =23
         numtets(27) =23
         numtets(28) =23
         numtets(29) =23
         numtets(30) =23
         numtets(31) =23
         numtets(32) =22
         numtets(33) =21
         numtets(34) =18
         numtets(35) =15
         numtets(36) = 1
      endif
      if(nlams.eq.37) then
         numtets( 1) = 5
         numtets( 2) = 7
         numtets( 3) = 8
         numtets( 4) =10
         numtets( 5) =11
         numtets( 6) =12
         numtets( 7) =13
         numtets( 8) =14
         numtets( 9) =15
         numtets(10) =16
         numtets(11) =17
         numtets(12) =17
         numtets(13) =18
         numtets(14) =19
         numtets(15) =19
         numtets(16) =20
         numtets(17) =21
         numtets(18) =21
         numtets(19) =22
         numtets(20) =22
         numtets(21) =22
         numtets(22) =23
         numtets(23) =23
         numtets(24) =23
         numtets(25) =24
         numtets(26) =24
         numtets(27) =24
         numtets(28) =24
         numtets(29) =24
         numtets(30) =24
         numtets(31) =24
         numtets(32) =24
         numtets(33) =22
         numtets(34) =22
         numtets(35) =21
         numtets(36) =14
         numtets(37) = 1
      endif
      if(nlams.eq.38) then
         numtets( 1) = 5
         numtets( 2) = 7
         numtets( 3) = 9
         numtets( 4) =10
         numtets( 5) =11
         numtets( 6) =12
         numtets( 7) =13
         numtets( 8) =14
         numtets( 9) =15
         numtets(10) =16
         numtets(11) =17
         numtets(12) =18
         numtets(13) =18
         numtets(14) =19
         numtets(15) =20
         numtets(16) =20
         numtets(17) =21
         numtets(18) =22
         numtets(19) =22
         numtets(20) =22
         numtets(21) =23
         numtets(22) =23
         numtets(23) =24
         numtets(24) =24
         numtets(25) =24
         numtets(26) =24
         numtets(27) =25
         numtets(28) =25
         numtets(29) =25
         numtets(30) =25
         numtets(31) =25
         numtets(32) =25
         numtets(33) =24
         numtets(34) =24
         numtets(35) =22
         numtets(36) =21
         numtets(37) =15
         numtets(38) = 1
      endif
      if(nlams.eq.39) then
         numtets( 1) = 5
         numtets( 2) = 7
         numtets( 3) = 9
         numtets( 4) =10
         numtets( 5) =11
         numtets( 6) =12
         numtets( 7) =13
         numtets( 8) =14
         numtets( 9) =15
         numtets(10) =16
         numtets(11) =17
         numtets(12) =18
         numtets(13) =19
         numtets(14) =19
         numtets(15) =20
         numtets(16) =21
         numtets(17) =21
         numtets(18) =22
         numtets(19) =22
         numtets(20) =23
         numtets(21) =23
         numtets(22) =24
         numtets(23) =24
         numtets(24) =24
         numtets(25) =25
         numtets(26) =25
         numtets(27) =25
         numtets(28) =25
         numtets(29) =25
         numtets(30) =26
         numtets(31) =25
         numtets(32) =25
         numtets(33) =25
         numtets(34) =25
         numtets(35) =25
         numtets(36) =23
         numtets(37) =20
         numtets(38) =15
         numtets(39) = 1
      endif
      return
      end
